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ABSTRACT 

Observations  of  positively  charged  particles  that  are 
generated  on  or  near  satellite  surfaces  have  been  made  on 
several  spacecraft.  This  thesis  postulates  sputtering  of 
the  satellite  surface  due  to  ambient  ion  impact  as  the 
generating  mechanism.  Calculations  are  made  using  the 
S igmund-Thompson  sputtering  theory  to  determine  the  response 
at  the  particle  detectors.  These  calculations  indicate  that 
surface  sputtering  creates  a  sufficient  flux  to  account  for 
the  observed  phenomena.  The  NASA  Charging  Analyzer  Program 
was  run  to  determine  the  trajectories  for  actually  observed 
particles.  The  calculated  trajectories  were  determined  to 
lead  to  the  spacecraft  surface,  again  indicating  that 
surface  emission  was  the  source.  The  sputtering  flux  as 
calculated  was  insufficient  to  cause  any  significant  short- 
term  damage  to  the  spacecraft,   beyond  thin  coating  erosion. 
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I .   INTRODUCTION 

A.   OVERVIEW 

Spacecraft  charging  is  the  development  of  a  potential 
difference  between  a  space  vehicle  and  its  plasma 
environment.  It  is  of  significance  for  scientific  purposes, 
since  the  potential  may  affect  measurement  of  the  ambient 
properties  of  the  environment,  and  for  practical  purposes, 
since  it  may  lead  to  anomalous  command  signals  and  physical 
damage  to  the  affected  spacecraft.  Although  predictions  of 
spacecraft  charging  had  been  made  previously,  the  first 
reports  of  high  level  spacecraft  charging  were  those  for 
Applications  Technology  Satellite  (ATS)  5  {Ref.  1}.  Similar 
events  were  subsequently  noted  on  ATS-6  and  P78-2  (SCATHA) 
by  various  observers  {Ref. 2}. 

Generally,  negative  charging  events  are  identified  by  an 
intense  peak  in  the  measured  ion  spectra  at  a  specific 
energy,  with  an  absence  of  ions  of  lesser  energies.  This 
intense  peak  is  most  clearly  seen  on  the  plot  of  the 
distribution  function  versus  energy,  as  will  be  displayed 
later  in  this  work.  This  peak  is  known  as  the  charging 
peak.  Such  observations  result  from  the  acceleration  of  the 
ambient  plasma  ions  into  the  detector  by  the  potential  drop 
between  the  plasma  and  the  spacecraft.  Occasionally,  ion 
fluxes   are  seen  at  energies  below  the   peak   energy,   which 


would  appear  to  violate  conservation  of  energy  for  a 
col 1  is ionless  plasma.  These  ion  distributions  have  been 
reported  previously,  and  tentatively  ascribed  to  sputtering 
from  spacecraft  surfaces  {Ref.  3}.  They  have  been  termed 
"spacecraft  generated  ions"  {Ref.  2}. 

Analysis  of  this  particular  phenomenon  is  of  practical 
interest,  as  these  ions  may  contaminate  the  charging  peak, 
leading  to  an  underestimation  of  the  satellite  potential. 
The  contamination  of  the  charging  peak  may  be  particularly 
important  active  experiments  in  low  earth  orbit, such  as 
SPACELAB .  Additionally,  since  we  believe  the  source  to  be 
sputtering,  the  flux  of  these  particles  is  an  indication  of 
the  damage  to  the  satellite  surface  caused  by  the  vehicle's 
environment.  Also,  the  life  expectancy  of  precision 
satellite  surfaces,  e.g.  optical  surfaces,  and  surface 
coatings  may  be  directly  affected  by  the  sputtering  rate 
appropriate  to  their  environment. 

B.   SCATHA  SPACECRAFT 

Most  of  the  experimental  data  and  calculations  performed 
in   this  thesis  are  based  on  the  P78-2  (SCATHA)   spacecraft. 

The  construction,   environment,   and  sensors  of  this  vehicle 
will  be  discussed  in  the  next  several  sections. 
1 .   Spacecraft  Construction 

The  U.S.  Air  Force  P78-2  spacecraft  was  launched  in 
January  1979  to  collect  data  for  a  study  of  Spacecraft 
Charging    AT    High   Altitudes    (SCATHA),    in    a    joint 


NASA/Department  of  Defense  Program.  The  satellite  body  was 
cylindrical  in  shape  with  a  length  and  diameter  of 
approximately  1.75  meters.  Seven  booms  were  deployed  in 
orbit  from  the  spacecraft  to  provide  isolation  for 
experiments  from  charging  effects  on  the  satellite  surface. 
Figure  1  depicts  the  SCATHA  spacecraft.  The  top  of  the 
satellite  in  the  figure  is  the  forward  end,  and  the  bottom 
is  aft.  The  experiments  and  communications  equipment  for 
the  satellite  are  primarily  located  on  equipment  decks  in 
the  center  of  the  cylinder.  The  University  of  California, 
San  Diego,  charged  particle  detector  was  on  the  forward  end, 
and  the  ion  gun  was  on  the  aft  end;  both  aro  described 
below.  The  outer  cylinder  surface  was  divided  into  three 
general  areas;  two  solar  arrays,  one  forward  and  one  aft, 
and  a  bellyband  between  them  to  provide  access  panels. 
Additionally,  sections  of  the  surface  were  covered  in 
various  materials,  such  as  kevlar  and  mylar,  for  use  in  the 
experiments.  (Ref.  4} 

The  spacecraft  was  constructed  of  special 
lightweight  materials,  as  its  orbit  was  at  the  maximum 
capability  of  the  launch  vehicle.  The  central  tube  of  the 
spacecraft  was  magnesium  and  the  equipment  decks  were 
aluminum.  The  solar  arrays  were  aluminum  core  with  a 
fiberglass  (Si02)  outer  face.  The  bellyband  panels  were 
covered   with  thermal  paint  and  second  surface   mirrors,  and 


Figure  1.  SCATHA  Satellite  {Ref.  4} 


acted  as  waste  heat  radiators.   The  forward  cylinder  end  was 
coated  with  gold,  except  for  some  sample  patches.   {Ref.  4} 
2 .   Spacecraft  Orbit  and  Environment 

The  vehicle  was  inserted  into  a  near-geosynchronous, 
elliptical  orbit  with  a  perigee  of  5.5  R«  and  an  apogee  of 
7.7  R« .  The  orbit  period  was  23.5  hours,  and  the 
inclination  to  the  equatorial  plane  was  7.8  degrees.  The 
satellite  axis  was  perpendicular  to  the  earth-sun  line,  and 
the  vehicle  had  a  spin  period  of  59  seconds.  The  satellite 
encountered  40  day  eclipse  periods  in  both  the  fall  and  the 
spring.  During  these  times,  the  spacecraft  was  shadowed  for 
approximately  1  hour  per  day.  Figure  2  gives  a  schematic 
of  the  orbit.  The  view  is  from  above  the  earth's  north 
pole .  {Ref . 5} 

The  environment  of  the  SCATHA  spacecraft  consisted 
of  two  main  regimes,  the  plasmasphere  and  the  plasmasheet. 
The  nature  of  these  regimes,  particularly  the  ion 
composition,  is  important  in  determining  the  magnitude  of 
the  vehicle's  interaction  with  the  environment,  for  example, 
the  sputtering  yield.  The  boundary  between  these  regimes 
(the  plasmapause)  is  well  defined,  but  varies  with  solar 
activity  and  solar  wind  parameters.  The  spatial  regions  are 
depicted  in  Figure  3.  The  view  depicts  the  northern 
hemisphere.  {Ref.  6} 

The  innermost  region  is  the  plasmasphere.  The 
density   is    generally   greater   than   10®   m~^.    and   the 


Figure  2.  SCATHA  Orbit  {Ref.4} 


Figure  3.  Regions  of  Space  near  Earth  {Ref.  6} 


temperature  is  less  than  1  eV,  The  dominant  ion  is  H"^,  with 
0"^  and  He"^  contributing  approximately  20°6  of  the  energy 
density.  The  source  of  these  ions  is  currently  considered 
to  be  the  ionosphere.  {Refs.  6,7,8} 

Overlapping  the  plasmasphere  is  the  ring  current. 
Tt  consists  of  a  hot,  thin  plasma,  with  densities  normally 
less  than  10®  m~^  and  temperatures  greater  than  10  keV. 
Again,  H"^  is  the  dominant  species  (93%  of  the  energy 
density)  while  0*  and  He*  are  typically  minor  constituents. 
During  quiet  times  the  principal  source  of  these  ions 
appears  to  be  the  solar  wind.  During  active  times,  a  large 
0*  constituent  indicates  ionospheric  sources.  {Refs. 
6,7,3} 

The  plasmasheet  is  the  region  outside  the 
plasmapause,  and  like  the  ring  current,  consists  of  a  hot, 
diffuse  plasma.  Densities  on  the  order  of  10*  -  10^  m~^  and 
temperatures  near  10  keV  are  common.  At  quiet  times,  the 
composition  is  90%  H"^  and  10%  0"^  and  He"^.  At  active  times, 
large  0"^  ratios  (>50%)  are  observed.  This  is  a  mixing  zone, 
containing  particles  of  both  terrestrial  and  solar  origin. 
Large  negative  potentials  are  found  on  spacecraft  in  the 
plasmasheet,  particularly  when  the  satellite  is  in  eclipse. 
Table  1  shows  the  various  ions  in  the  magnetosphere,  and 
their  sources  .  {Refs  .  6,7,8} 

An  interesting  feature  of  these  areas  is  the 
variability   of  the  ion  composition  of  the  plasmas   observed 


during  magnetically  active  periods.  In  the  ring  current,  for 
example,  it  has  been  observed  that  0"^  contributes  a  minimum 
of  23%  of  the  storm  time  energy  density,  and  may  contribute 
much  more.  Also,  the  He*  contribution  increases,  and 
adiabatic  acceleration  of  these  solar  wind  particles  may 
leave  them  with  energies  up  to  32  keV.  Table  2  gives  some 
typical  values  for  storm  time  ring  current  composition. 
{Refs.   7,8} 

The  plasmasheet  region  has  a  direct  source  of  0*,  in 
field  aligned  beams  streaming  from  the  ionosphere.  The 
H'^/0"^  ratio  of  the  beams  has  been  found  to  be  6/1  [Eef. 
8}.  Additionally,  it  has  been  observed  that  in  the  L  =  G  - 
3  region  (e.g.,  the  SCATHA  orbit),  during  periods  of  high 
magnetic  activity  (Kp~5),  the  0"^  density  becomes  comparable 
to  that  of  H"^.  (The  Mcllwain  'L'  parameter  is  defined  by 
R=L  cos^S,  where  R  is  the  distance  in  earth  radii,  and  9  is 
the  magnetic  latitude.)  K^  is  a  measure  of  the  general 
level  of  magnetic  activity  caused  by  the  solar  wind.  The  Kp 
scale  is  quas i -logar ithmic,  and  runs  from  0-9,  with  0  the 
low  end  and  9  the  high.  On  several  occasions,  the  0"^ 
density  has  been  found  to  reach  70%.  Table  3  lists  some 
values  for  the  storm  time  magnetospher ic  bulk  plasma 
composition.  {Refs.  5,7,0} 
3  .   Instruments 

The  data  for  this  experiment  came  primarily  from  the 
University   of   California,    San   Diego   auroral   particles 


experiment,  SC-9 .  This  experiment  was  composed  of  five 
detectors;  two  pairs  of  rotating  ion  and  electron  detectors, 
and  one  fixed  ion  detector.  Figure  4  depicts  the  detector 
arrangement.  The  rotating  detectors  scanned  in  orthogonal 
planes  and  were  designated  the  Hi  set  and  Lo  set.  The  Lo 
detector  rotated  from  -20  deg  to  200  deg,  with  70  deg  along 
the  spacecraft  spin  axis.  The  Hi  detector  scanned  the  same 
range,  but  it  was  more  symmetrical,  with  its  midpoint  on  a 
line  parallel  to  the  spin  axis.  Thus,  it  could  not  depress 
more  than  20  degrees  below  the  forward  end  plane.  The  fixed 
detector  looked  radially  away  from  the  spacecraft , in  the 
same  plane  as  the  Lo  detector.  The  detectors  had  an  angular 
resolution  of  5  by  7  deg.  The  geometry  is  illustrated  in 
Figure  5.  {Ref.  4} 

The  Hi  set  covered  an  energy  spectrum  from  1  eV  to 
81  keV.  The  Lo  rotating  set,  and  the  fixed  detector, 
scanned  from  1  eV  to  2  keV.  For  all  detectors,  a  complete 
energy  scan  required  16  seconds  and  was  covered  in  a  series 
of  64  logarithmic  steps.  The  energy  resolution  at  each  step 
was  approximately  20*%.  Additionally,  the  detectors  could  be 
ordered  to  dwell  at  fixed  positions  and/or   energies . {Ref . 4 } 

The  detectors  were  composed  of  three  parts;  an 
electrostatic  analyzer,  an  electrostatic  focussing  lens,  and 
a  spiraltron  particle  sensor.  The  analyzer  provided  the 
energy  differentiation  for  the  system  through  curved  plates 
and  an  applied  voltage.    The  lens  focussed  the  particles  on 
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the  sensor  center  by  means  of  two  grids,   one  at  ground   and 
the  other  at  the  analyzer  potential.    The  sensor  was   rated 
to  10"^  counts  per  second.  {Ref.  4} 
4 .   Ion  Gun 

An  ion  gun  was  installed  on  the  SCATHA  spacecraft  to 
investigate  the  efficiency  of  an  ion  emission  system  in 
modifying  satellite  potentials.  In  particular,  it  could  be 
used  to  develop  negative  voltages  on  the  vehicle.  The 
experiment  utilized  xenon  gas  ionized  by  cathode  discharge 
and  accelerated  by  either  a  1  or  2  kV  potential  drop.  The 
beam  current  could  be  varied  incrementally  from  0.3  mA  to 
2.0  mA .  The  package  also  included  an  electron  source  that 
could  be  configured  to  neutralize  the  particle  beam  or  to 
provide  an  electron  beam.  (Ref.  9} 

C.   OTHER  SATELLITES 

Although  SCATHA  was  the  primary  source  of  data  for  this 
thesis,  two  other  satellites  provided  significant 
information.  These  satellites,  ATS-6  and  ISEE-1,  will  be 
briefly  discussed  in  the  following  paragraphs. 

1.   ATS-6 

ATS-6  was  a  large  three-axis  stabilized  spacecraft 
placed  in  geosynchronous  orbit  in  1974.  Its  major  features 
are  shown  in  Figure  6  .  The  large  antenna  was  significant, 
since  it  fostered  differential  charging  by  shadowing  other 
spacecraft  surfaces.  The  satellite  had  an  environmental 
measurements   package   above   the   spacecraft   dish.     This 
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package  contained  rotating  particle  detector  sets  similar  to 
those  found  on  SCATHA.  {Ref.2} 
2.   I  SEE 

The  ISEE-1  satellite  was  launched  in  1977  as  part  of 
a  three  spacecraft  mission  to  investigate  the  magnetospher ic 
plasma.  With  this  in  mind,  the  satellite  was  designed  to 
minimize  absolute  and  differential  charging  by  using 
conducting  materials  and  low  impedance  connections  to 
ground.  Regardless,  the  satellite  was  observed  to  charge  on 
several  occasions.  The  spacecraft  was  equipped  to  measure 
electric  fields,  plasma  waves,  and  plasma  composition.  The 
most  important  of  these,  for  our  purposes,  is  the  plasma 
composition  experiment  which  contained  an  ion  mass 
spectrometer.  {Ref.  10} 

D.   SPACECRAFT  CHARGING 

The  theory  of  spacecraft  charging  will  now  be  reviewed, 
as  differential  charging  will  play  a  significant  role  in 
later  discussions.  In  general,  a  probe  immersed  in  a  plasma 
will  develop  some  potential  relative  to  that  medium.  This 
"floating"  potential  is  determined  by  the  balance  of  the 
currents  incident  upon  the  probe.  The  primary  currents  in 
deciding  the  potential  are;  the  ambient  plasma  electron  and 
ion  flux,  the  photoelectron  emission  due  to  sunlight,  and 
the  fluxes  due  to  secondary  and  backscattered  electrons. 
These  currents  are  schematically  illustrated  in  Figure  7. 
The   left   half   of   the   figure   depicts   typical    eclipse 
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conditions   (negative   satellite),   and  the  right   half   the 
currents  present  in  sunlight.  {Refs.  11,12,13} 

In  daylight,  the  dominant  current  is  that  of  the 
photoe lectrons  which,  although  material  dependent,  has 
typical  saturation  current  densities  in  the  range  5-50 
pA/m2.  This  current  density  overshadows  that  of  the  ambient 
plasma,  which  is  generally  below  1  \lA/m^  for  electrons,  and 
.1  ]lA/m-  for  ions.  The  value  of  the  secondary  electron 
emission  current  is  highly  dependent  on  both  the  material 
and  the  electron  temperature.  It  is  this  material 
dependence  which  allows  differential  charging,  even  in 
eclipse.  A  normal  secondary  emission  curve  will  have  a 
yield  greater  than  one  in  some  middle  energy  range 
(typically  50-500  eV),  and  yields  less  than  one  at  higher 
and  lower  energies.  The  point  in  the  high  energy  range 
where  the  yield  decreases  to  one  is  called  the  cross-over 
point  and  is  critical  in  determining  the  effect  of  the 
secondary  emission  current.  The  energy  at  which  this  point 
occurs,  and  the  peak  yield,  varies  with  target  composition. 
Low  energy  ambient  distributions  (T<5  keV}  impacting  the 
spacecraft  will  create  large  secondary  fluxes,  in  fact 
greater  than  the  ambient  electron  flux,  and  the  spacecraft 
will  charge  positively.  High  energy  electron  distributions 
(T>5  keV)  will  produce  lower  secondary  emission  yields,  and 
the  vehicle  will  gain  negative  charge.  Thus,  the  potential 
of   the  satellite  is  determined  by  the  high  energy   electron 
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flux.  The  threshold  electron  temperature  for  which  charging 
may  occur  is  typically  in  the  5-10  keV  region.  There  is  an 
additional  citeria,  the  critical  energy,  which  also  must  be 
met.  This  is  an  upper  energy  bound  on  the  distribution 
function  resulting  from  magnetospher ic  convection  processes. 
It  is  also  termed  the  Alfven  boundary,  and  must  exceed  15-20 
keV  for  charging  to  occur.  {Refs.  11,12,13,14,15,16} 

A  simplistic  equation  to  represent  the  situation  may  be 
given  by: 

(1-1)  Iri«-t:  —       Iamb   +   I  » •  <=   +   I  p  hi  o  t  o 

Thus,  during  daylight,  when  the  photoelectron  current 
dominates,  the  net  current  is  positive  (toward  the 
spacecraft),  and  the  spacecraft  will  start  to  charge 
positively.  As  this  happens,  the  flow  of  ions  to,  and 
electrons  from,  the  spacecraft  will  be  inhibited,  and  the 
value  of  In.t  will  decrease.  When  the  value  reaches  zero, 
the  currents  are  balanced  and  the  spacecraft  will  remain  at 
that  potential  until  the  balance  is  disturbed.  For  eclipse 
events,  Ip^oco  is  zero,  and  the  balance  is  effectively 
between  the  ambient  and  secondary  electrons.  For  a  low 
energy  electron  distribution,  the  potential  is  restricted  to 
low  values  (+2  to  +5  V),  while  for  high  energies,  the 
potential  may  reach  -1  to  -20  kV.fRefs.  2,11} 

In  addition  to  absolute  charging,  a  spacecraft  may 
charge  differentially.  That  is,  separate  insulated  portions 
of    the    spacecraft   surface   may   charge    to    different 
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potentials.  If  two  insulators  are  exposed  to  the  same 
environment,  but  have  different  secondary  emission 
characteristics,  then  their  ability  to  shed  the  excess 
charge  in  the  form  of  a  current  will  differ,  and  a  potential 
will  develop  between  them.  Also,  if  two  portions  are 
exposed  to  differing  environments,  they  will  each  charge  to 
different  potentials.  The  environment  may  differ  over  the 
range  of  a  spacecraft  in  a  number  of  ways.  Frequently, 
pitch  angle  anisotropies  in  the  ambient  plasma  may  cause 
different  fluxes  to  impact  separate  spacecraft  regions. 
Additionally,  if  one  portion  of  the  spacecraft  shadows 
another  continually  (see  ATS-6),  the  shadowed  area  will  not 
emit  photoelectrons  and  will  charge  negatively  compared  to 
the  rest  of  the  satellite.  Also,  in  lower  orbits,  an 
electric  field  induced  in  the  satellite  frame  of  reference 
(E=vxB)  may  cause  variations  in  the  ambient  particle 
distributions.  {Refs.  2,12,13} 

Differential  charging  is  believed  to  be  the  cause  of 
numerous  control  logic  upsets  on  various  spacecraft.  As  the 
differential  potential  increases,  there  will  be  a  tendency 
for  it  to  discharge  in  the  form  of  an  arc.  This  arc  creates 
an  electromagnetic  pulse  which  is  read  by  the  control 
circuits  as  a  command  signal.  For  large  discharges  there 
may  even  be  physical  damage  or  destruction  to  the  arcing 
areas.  (Refs.  17,18} 
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Spacecraft  charging  may  be  controlled  by  active  and 
passive  measures.  Active  trials  using  ion  and  electron  guns 
were  carried  out  on  ATS-S  and  SCATHA  with  some  success. 
Passive  attempts,  such  as  those  on  the  ISEE  satellites, 
using  low  impedance  connections  and  conducting  surfaces  were 
also  moderately  successful.  {Refs.  10,19,20} 

E.   SPUTTERING 

The  general  theory  of  sputtering  is  very  significant  to 
this  thesis,  as  it  is  postulated  to  be  the  source  of  the 
observed  spacecraft  generated  ions.  Sputtering  is  defined 
as  the  ejection  of  material  from  solid  surfaces  under  ion 
bombardment.  The  observed  features  of  the  sputtering  yield 
(number  of  target  atoms  sputtered  per  incident  ion)  are;  a) 
its  dependence  on  atomic  number,  b)  its  correlation  with  the 
periodic  table,  c)  its  correlation  with  the  heat  of 
sublimation  of  solids,  d)  its  relation  to  the  efficiency  of 
momentum  transfer,  e)  and  its  dependence  on  incident  ion 
energy,  such  that  it  increased  to  a  maximum  and  then 
decreased.  Figure  8  is  a  plot  of  the  sputtering  yield  for 
aluminum  as  a  function  of  incident  ion  energy  and  mass  {Ref. 
21} .  We  see  that  the  yield  maximum  increases  with  ion  mass, 
and  shifts  to  higher  energies.  These  observations  led  to 
theories  based  on  atomic  collision  cascades.  Such  cascades 
are  created  by  one  incident  ion  impinging  on  the  lattice  and 
transferring  energy  to  other  atoms,  which  in  turn 
participate  in  collisions.    These  theories  are  based  on  the 
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Figure  8..  Sputtering  Yield  for  Aluminum  {Ref.  21} 
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simple  elastic  collisions  of  classical  mechanics,  and 
therefore  satisfy  the  observed  dependence  on  atomic  number 
and  momentum  transfer.  The  correlation  to  the  heat  of 
sublimation  is  explained  by  assuming  the  effect  of  some 
surface  binding  energy.  The  maximum  in  the  energy 
dependence  occurs  by  postulating  some  optimum  energy 
deposition  depth,  beyond  which  the  ability  of  the  cascade 
particles  to  reach  the  surface  decreases.  {Refs.  21,22,23} 

The  most  comprehensive,  detailed,  and  successful  of  the 
collision  cascade  theories  is  the  Sigmund  -  Thompson  theory. 
This  theory  is  actually  two  theories  developed  separately, 
but  they  complement  one  another  by  explaining  different 
facets  of  sputtering.  The  portion  developed  by  Sigmund 
treats  the  sputtering  yield  of  an  amorphous  target. 
Thompson's  theory  explains  the  energy  distribution  of  the 
sputtered  atoms.  The  development  of  both  these  theories 
will  be  briefly  discussed,  followed  by  a  comparison  with 
experimental  results.  {Refs.  22,23,24,25} 

1 .   The  Sigmund  Theory 

The  Sigmund  theory  focusses  on  the  amount  of  energy 
deposited  in  the  surface  layers  of  the  solid  as  the  driving 
factor  for  the  yield.  The  basic  expression  for  the 
sputtering  yield  is  then  S=rF,  where  r  is  a  constant 
depending  only  on  the  properties  and  state  of  the  target, 
and  F  is  a  function  describing  the  interaction  of  the  two 
particles.    In  order  to  evaluate  this   expression,   Sigmund 
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chose  as  his  starting  point  the  sputtering  of  an  amorphous 
target,  and  used  transport  theory  to  describe  the  collision 
cascade   in  random  media.   {Refs.   23,24} 

The  theory  assumes  a  planar  surface  and  an  atom 
starting  its  motion  at  time  t  =  0  and  position  x  =  0.  The 
initial  quantity  of  interest  is  G  (  x,a/i3  ,  v,  t  )d^Va  dx  .  This 
is  the  average  number  of  atoms  moving  at  time  t  in  a  layer 
(x,dx)  with  a  velocity  (v's/d^Va)  due  to  an  impact  by  a 
single  atom  at  velocity  v.  This  quantity  is  then  used  in  a 
Boltzmann  transport  equation  that  equates  the  initial  number 
G  with  a  later  G  which  is  dependent  on  the  probability  of 
collision.  The  equation  is  then  integrated  to  eliminate 
some  of  the  variables.  Two  functions  are  defined  to 
simplify    the    Integration. 


(1-2) 


(1-3) 


F(x,V(a,V)      =      r  G(X,Va,^,t)dt 
H(X,V)   =   P  F( X, Va, V) I Vax I d^Va 


J 


The  function  F  is  the  total  number  of  atoms  which  penetrate 
the  plane  x  with  a  velocity  (v'a^d^va)  during  the  collision 
cascade  development.  H(x,v)  is  the  sputtering  yield  at  the 
plane  x  for  a  source  of  sputtered  particles  at  x  =  0.  This 
function  may  be  defined  for  either  forward  or  backward 
sputtering.  For  our  purposes,  we  will  choose  backward 
sputtering.  Backward  sputtering  is  emission  directed 
opposite  to  the  initial  direction  of  inc idence .  { Re f s  .  23,24] 
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The  remaining  expression  in  H  is  expanded  in  terms 
of  Legendre  polynomials  and  then  transformed  into  a  moment 
equation  by  multiplying  each  term  by  x"  and  integrating  over 
X.  With  the  assumption  that  the  electronic  losses  to 
electrons  are  small  and  isolated  from  nuclear  collisions, 
the  electronic  and  collisional  cross  sections  may  be 
separated  {Ref.26}.  This  yields  a  final  expression  which 
includes  the  electronic  stopping  cross  section,  and  the 
differential  cross  section  for  elastic  scattering.  {Refs. 
23,24} 

The  next  step  is  to  determine  the  functions 
describing  the  electronic  and  elastic  cross  sections.  The 
electronic  cross  section  used  is  S«  =  K  E"*,  where  the 
constant  K  is  determined  using  Thomas-Fermi  arguments 
{Ref.27}.  The  elastic  cross  section  is  derived  from 
classical  mechanics,  and  a  power  approximation  of  the 
Thomas-Fermi  interatomic  potential.  The  equation  is  then 
given  by; 


(1-4) 


dCT(T)   =   C   E-m-p-X-m   (3^1 


where  C  is  a  constant,  T  is  the  transferred  energy,  and  E 
the  initial  energy.  The  variable  m  comes  from  the  power 
approximation.  m=l  corresponds  to  Rutherford  scattering,  and 
m=0  approximates  scattering  from  a  Born-Mayer  potential.  If 
these  values  are  substituted  in  the  expression  from  the 
previous  paragraph,  the  following  is  achieved; 


(1-5)      H(x,E,n)  =  3_  F(x.E,B) 

4ii=  N  Cb  Ub 


where  3  is  the  direction  cosine  of  the  ejection  vector,  and 
F(x,E,f3)  is  a  function  describing  the  distribution  of  the 
energy  deposited  in  the  solid.  Ca  is  a  constant,  and  Ua  is 
the  binding  energy.  In  general,  the  function  F(x,E,3)  is 
found  to  obey  the  following  relation; 


(1-6 


F(x,E,0)  =  a  N  Sr^(E) 


where  a  is  a  dimens ionless  constant  dependent  on  the  ratio 
M2/M1,  N  is  the  atom  density,  and  Sr»(E)  is  the  nuclear 
stopping  power,  given  by  Sn(E)  =  P  T  da.  Additionally,  we 
now  restrict  ourselves  to  sputtering  in  the  unit  normal 
direction.  The  function  a  is  illustrated  in  Figure 
9  .  {Refs.23,24} 

Using  all  of  the  previous  relations,  we  can  then 
define  some  general  expressions  for  the  sputtering  yield. 
For  low  energies  (E<lkeV),  we  have 


(1-7)      S(E,3)  =  3_  a   4M^M^    E_ 

4Tt2    (Mi  +  Ms)^  Ua 


Sn(E)  in  this  equation  was  evaluated  by  integrating  dcr  with 
a  value  of  m=0.  The  Born-Mayer  potential  is  more  accurate 
at  lower  energies  than  the  Thomas-Fermi  approximation.  The 
expression  for  higher  energies  is; 


(1-8) 


S(E,(3)  =  3  .  56  a 


'  +  Z2^^3  ) 
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Figure    9.    a(Mi/Mi)    vs.    Mi/Mi    {Ref.    23} 


25 


The  function  Sr»(E)  is  the  reduced  stopping  power.  This  is  a 
universal  function  for  Thomas-Fermi  interactions  {Refs. 
24,27}.   The   variable  E    is   the  reduced  energy,   defined  by 


e  = 


aM^E 


,  where   Zi  and    Zs    are   the   atomic 

(ZiZ^e^ (MX+M2) ) 
number,   and  a  is  the  Lindhard  screening  radius.    Tabulated 

values  for  Sn(£)  are  shown  in  Table  4.   The  sputtering  yield 

dependence   on   0  was  roughly  determined  to   be   (cos(9))~^. 

This    dependence   was   derived   empirically    by    Sigmund. 

{Refs. 22, 23, 24, 27} 

2 .   The  Thompson  Theory 

The   Thompson   Theory   was   a   complete   theory    of 

sputtering   published  in  1968.    The  portion   involving   the 

determination   of  the  sputtering  yield  was   overshadowed   by 

the   Sigmund  Theory,   but  the  energy  spectrum   determination 

has  endured.   Thompson  started  his  derivation  by  assuming  an 

infinite   solid  with  recoils  at  energy  E2  from  some   primary 

event.    Then,   with  a  density  of  recoil   atoms  q(Ez)dEi  per 

unit  time,   each  generating  their  own  collision  cascade,  the 

total  number  of  atoms  slowing  down  through  some  energy  range 

E'  is 

(1-9)      n(E')  =  P-  q(E=)  M(E2,E')dE= 

J  h;  ' 

The  mean  rate  of  energy  loss  is  given  by  vdE'/dx,   where  v  is 
the  velocity.   Assuming  an  isotropic  distribution,   one   may 
def  ine 


(1-10) 


r 


n(E' ,r  )dE'd2'  =  p-  aiEz,)U(E^  .E'  )    dE^  dE  '  d^' 

V  dE'/dx  4Tt: 
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From  this,  the  flux  may  be  determined. 


(1-11)     *(E' ,r  )dE'dii'=  v  n  (  E  '  ,r  )  cosGdE  '  dii ' 


The  theory  now  assumes  that  the  infinite  solid  is 
cut  in  half  and  the  flux  is  observed.  Also,  some 
expressions  are  needed  for  the  functions  above.  H(E2",E')  is 
the  number  of  displaced  atoms  and  may  be  approximated  by 
(BEz/E',  where  (3  is  a  constant  of  order  unity.  Similarly, 
dE'/dx  may  be  replaced  with  E'/D,  where  D  is  the  interatomic 
spacing.   Using  these  relations  we  obtain 
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(1-12)     ^(E' ,r  )dE'dil'=i3_D  p-q(E=)E2dE2  cos9  dfl'  dE ' 

4  Ti: 


This  function  <!>  is  the  flux  inside  the  surface,   and  must  be 

transformed  to  that  observed  outside.    This  is  accomplished 

by  assuming  a  binding  force  normal  to  the  surface  which  will 

decrease  the  normal  kinetic  energy,   but  leave  the   parallel 

component  unaffected.   Thus,  the  effect  will  be  a  bending  of 

the  particles  trajectory.    With  the  energy   equations,   and 

the  relationship  between  the  inner  and  outer  angles  in  hand, 

we  come  to  the  final  expression 

(1-13)     'i'(E.0)dEd2=  B  D  cos0     P-  q(E2)E2dE2  dEd^l 

4ti(  l+Ua/E)  ^E^Jb- 

In  this  equation,  Ua  is  the  binding  energy,  E  is  the  emitted 
particle  energy,  and  0  is  the  angle  of  the  trajectory  from 
the  surface  normal.  The  chief  features  are  the 
proportionality   to   E/(E+U0)^   and   the   cos*    dependence. 
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Figure  10  is  a  plot  of  this  function  versus  energy.  For 
high  energies,  the  spectrum  will  fall  off  as  l/E^,  and  for 
low  energies,  the  binding  energy  will  have  an  impact.  The 
maximum  of  this  curve  occurs  at  U0/2.  Figure  11  is  a  plot 
of  the  measured  energy  d istr ibut ion, with  the  model  overlaid. 
It  can  be  seen  that  up  to  approximately  1000  eV,  the  model 
matches  the  observed  spectrum  quite  closely.  Since  the 
majority  of  the  sputtered  particles  are  in  the  energy  region 
that  follows  the  1/E^  curve,  this  theory  is  quite  effective 
in  describing  the  data.  {Refs.  22,24} 

3 .   Comparison  of  Theory  with  Experiment 

To  evaluate  the  applicability  of  these  theories,  we 
must  determine  how  closely  they  recreate  experimental 
results.  Therefore,  we  will  look  at  the  sputtering  yield  as 
a  function  of  ion  energy,  angle  of  incidence,  incident  and 
target  masses,  and  surface  binding  energy.  The  energy 
dependence  of  the  sputtering  yield  is  contained  in  the 
energy  dependence  of  the  stopping  power.  Calculated  yields 
using  equations  1-7  and  1-8  are  compared  with  experimental 
values  in  Figure  12  for  copper  {Ref.  23}.  The  solid  line 
is  the  high  energy  result,  and  the  dotted  line  is  the  plot 
for  the  low  energy  equation.  The  experimental  results  are 
marked  by  the  various  geometric  shapes.  In  addition, 
calculated  and  measured  yields  for  various  materials  are 
listed  in  Table  5  {Ref.  23}.  It  can  be  seen  that  the 
correlation  is  reasonable,   normally  within  a  factor  of  two. 
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Figure  11.  Energy  Distribution  of  Sputtered  Gold  Atoms 
{Ref.  23} 
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Some  of  the  experimental  yields  are  questionable,  as  the 
sputtering  yield  is  known  to  have  a  dose  dependence.  The 
Sigmund  Theory  also  tends  to  overestimate  sputtering  at  very 
low  energies  ( <100  eV).  In  most  cases,  experimental  data 
indicates  that  the  yield  goes  to  zero  in  the  50-100  eV 
region . {Refs  .  22,23,24} 

The  sputtering  yield  is  a  function  of  Mi  directly 
and  in  the  factor  aCMs/Mi).  Due  to  the  above  mentioned  dose 
effects,  yield  data  for  masses  are  generally  normalized  to 
the  self  sputtering  yield  or  the  argon  sputtering  yield. 
Data  for  various  targets  are  presented  in  Figure  13  {Ref. 
23}.  The  top  plot  displays  results  for  impacts  on  silicon. 
The  mass  of  the  sputtering  ion  in  AMU  is  along  the 
horizontal  axis  and  the  yield  of  the  ions  normalized  to  the 
argon  sputtering  yield  on  the  vertical  axis.  The  other 
plots  are  similar  for  copper  and  gold.  The  agreement  with 
theory  for  silicon  is  quite  good,  but  decreases  for  the 
heavier  targets.  Sputtering  yields  for  H*  and  0"^  on  silicon 
are  illustrated  in  Figure  14.  Since  in  general  we  will  be 
discussing  silicon  or  silicon  dioxide  as  a  target,  we  will 
disregard  the  lack  of  agreement  for  heavier  targets  in  the 
model  which  we  develop.  {Refs.  22,23,24} 

The  angular  dependence  of  the  sputtering  yield  is 
considered  t.o  be  proportional  to  (l/cos8)".  The  exponent, 
n,  is  a  variable  quantity  ranging  from  a  little  less  than  1 
to   1.7.   For   the  situations  we  will   be   discussing,   with 
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Figure  12.  Calculated  and  Experimental  Sputtering  Yields 
for  Copper  {Ref.  23} 
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Figure    14.    Typical    Sputtering   Yields    for    H+   and    0+ 
on   a    Silicon   Target 
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M=/Mi  >  1,  the  value  of  n  is  generally  taken  to  be  1.7. 
Figure  15  is  a  plot  of  the  dependence  on  angular  incidence 
from  experimental  data  {Ref.  23}.  The  dependence  ranges  up 
to  an  incident  angle  of  70  degrees.  Beyond  this  point,  the 
sputtering  yield  decreases  rapidly  to  0.  If  we  attempt  to 
model  using  only  the  0-70  degree  range,  we  will  probably  be 
underestimating  the  yield  by  approximately  20^0.  {Refs. 
22,23,24} 

The  yield  depends  directly  on  the  binding  energy, 
thus  the  determination  of  this  energy  is  quite  important. 
Sigmund  suggested  that  the  appropriate  value  for  the  binding 
energy  would  be  the  sublimation  energy  of  the  material. 
This  generally  leads  to  sputtering  yield  values  that  are 
slightly  high.  Recent  work  has  given  theoretical  backing 
for  a  higher  value,  but  no  definite  numerical  quantity  was 
put  forward  {Ref.  28}.  It  has  been  suggested  that  the  value 
providing  the  best  results  is  17/10  of  the  sublimation 
energy  {Ref.  29}.  If  this  value  is  used  to  recalculate  the 
numbers  in  Table  5,  we  get  the  results  in  Table  6.  It  is 
obvious  that  in  general,  the  calculated  yields  in  Table  6 
are  much  closer  to  the  experimental  yields,  than  those 
calculated  for  Table  5.  For  calculations  in  this  thesis,  we 
will  use  the  larger  value  of  the  binding  energy.  (Refs. 
22,23,24} 
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Figure  15.  The  Dependence  of  the  Sputtering  Yield 
on  Angle  of  Incidence  {Ref.23} 
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4 .   Limitations  of  the  Sigmund  Theory 

The  theory  was  devised  for  amorphous  monatomic 
targets  bombarded  by  atomic  ions.  It  is  not  so  readily 
applicable  to  the  bombardment  of  molecular  solids,  as  will 
be  necessary  in  later  calculations.  Since  several  of  the 
reference  data  points  are  concerned  with  bombardment  with 
molecular  oxygen,  it  will  be  relevant  to  discuss  this  type 
of  impact.  The  case  of  molecular  ions  on  molecular  targets 
may  be  treated  through  computer  simulation,  but  it  is  an 
extensive,  complicated  project,  as  the  probabilities  for 
each  collision  combination  must  be  calculated.  In  most 
cases,  it  is  more  appropriate  to  utilize  experimental  data 
if  it  is  available.  Such  data  is  often  difficult  to 
obtain.  Only  recently  has  an  effort  begun  to  systematically 
tabulate  the  known  experimental  results  and  compare  them 
with  the  theory.  The  presently  tabulated  data  concerns 
atomic  ions  on  monatomic  targets  {Ref.  30}.  A  great  deal  of 
data  is  also  distributed  randomly  through  the  literature  on 
other  ion-target  combinations. 

An  additional  bit  of  information  is  available  for 
the  yields  of  monatomic  molecules.  It  has  been  found  for 
several  atoms  that  the  sputtering  yield  per  atom  for 
molecular  ions  is  slightly  higher  than  that  for  atomic  ions 
{Ref.  31}.  The  results  are  listed  in  Table  7.  The  yields 
are  greater  by  approximately  a  factor  of  two,  while  for 
silicon   they   range   from   1.15   to   1.30.     There   is    a 
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difference,  but  it  is  not  drastic.  While  it  is  difficult  to 
accurately  model  the  more  complicated  situations,  it  should 
be  possible  to  obtain  some  approximate  results  for  monatomic 
molecules  that  are  within  a  factor  of  two  of  the  actual 
results . 

F.   NASA  CHARGING  ANALYZER  PROGRAM 

The  NASA  Charging  Analyzer  Program  (NASCAP)  was 
developed  to  accurately  model  the  dynamics  of  spacecraft 
response  to  realistic  plasma  environments.  NASCAP  can 
effectively  simulate  the  charging,  in  both  laboratory  and 
magnetosphere  environments  of  objects  that  are 
geometrically,  materially,  and  electrically  complex.  The 
program  utilizes  a  timestep  procedure,  calculating  quasi- 
static  steps.  The  dynamics  are  driven  by  charge 
accumulation  from  external  sources,  charge  depletion,  and 
conduction  in  dielectrics.  Each  timestep  includes  a  full 
three  dimensional  electrostatic  potential  calculation.  The 
calculation  of  Poisson's  equation  over  a  grid  is  time- 
staggered  with  a  procedure  in  which  incident  charged 
particle  fluxes,  leakage  currents,  emission  currents,  and 
induced  spacecharge  effects  are  computed  based  on  the 
current  quasistatic  conditions.  NASCAP  is  limited  to 
situations  where  the  Debye  length  is  long  compared  to  the 
dimensions  of  the  examined  object.  This  condition  is  met 
for  the  SCATHA  spacecraft.   NASCAP  can  output  a  time  history 
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of  spacecraft  charging,  potential  contours,  charge  contours, 
current  contours,  and  particle  trajectories. 
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II .   OBSERVATIONS 

A.   PREVIOUS  OBSERVATIONS 

Data  taken  from  three   satellites,   AT3-5,   ATS-6,   and 
ISEE-1,    will   be   discussed,   prior   to   presenting   P78-2 
(SCATHA)   results.     Observations   of   the   three    primary 
phenomena  (triangle  peak,   shadow  peak,   diffuse  background) 
on  satellites  prior  to  SCATHA  will  be  illustrated. 

1.   ATS-5 

The  presence  of  ions  below  the  charging  peak  was 
first  reported  for  ATS-5  {Ref.l}.  ATS-5  data  taken  in 
eclipse  in  1969  and  1970  showed  ion  fluxes  in  a  broad  range 
of  energies  below  the  charging  peak.  A  typical  example  of 
such  data  is  shown  in  Figure  16.  The  data  is  presented  in 
grey  scale  spectrogram  format,  a  common  method  of  displaying 
the  data  concisely.  The  figure  is  vertically  separated  into 
three  components.  The  top  primarily  displays  magnetic  field 
information,  unnecessary  for  our  purposes.  The  remaining 
two  sections  indicate  the  count  rates  of  electrons  and  ions 
in  the  detector  versus  time  and  energy.  The  horizontal  axis 
defines  increasing  time,  and  in  this  case  is  labelled  with 
the  hours  of  day  274  of  1970.  The  vertical  axis  is 
associated  with  the  energy  of  the  particles.  The  minimum 
(3i50  eV)  of  the  energy  axis  is  located  between  the  two 
sections,   and   energy  increases  upward  for   electrons,   and 
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Figure  16.  ATS-5,  Spectrogram,  Day  274  of  1971 
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downward  for  ions  to  a  maximum  of  50  keV.  The  magnitude  of 
the  detected  countrates  is  indicated  by  the  shading  of  the 
figure  from  black  to  white,  with  black  corresponding  to  zero 
counts,  and  white  with  maximum.  If  the  detector  is 
saturated  at  a  certain  time  in  an  energy  channel,  the 
spectrogram  will  contain  a  black  spot  at  that  time  and 
energy.  A  feature  such  as  this  can  usually  be  distinguished 
from  a  zero  count  rate  because  saturation  points  are 
normally  surrounded  by  intense  count  regions,  revealed  as 
bright  white  areas. 

The  data  displayed  in  Figure  16  are  from  day  274  of 
1970.  The  satellite  enters  a  region  of  hot  plasma  at  04  5  0, 
as  indicated  by  the  increased  intensity  of  high  energy 
(>4500  eV)  electrons.  The  satellite  charging  peak  may  be 
seen  as  the  intense  white  region  in  the  proton  spectrum 
between  0620  and  0720.  From  the  energy  of  this  region,  we 
can  infer  a  satellite  potential  of  approximately  .75-2  kV. 
The  black  spots  on  the  peaks  are  an  example  of  how  detector 
saturation  is  displayed.  The  broad  spectrum  of  ions  is 
evident  in  the  energy  bands  below  the  charging  peak  as  the 
grey  area  that  clearly  extends  from  1  kV  to  about  300  V,  and 
may  even  continue  to  lesser  energies.  These  particles  did 
not  appear  during  every  charging  event,  so  the  source  was 
difficult  to  determine  and  no  complete  explanation  was  put 
forward . 
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2.   ATS-6 

ATS-6  also  showed  evidence  of  non-ambient  source 
contamination  of  the  ion  data.  Secondary  or  shadow  peaks 
were  observed  at  energies  below  the  inferred  satellite 
potential.  An  example  of  such  data  is  illustrated  in  Figure 
17.  This  figure  is  a  grey  scale  spectrogram  similar  to  the 
one  shown  in  Figure  15.  The  data  presented  are  from  the 
North-South  detector  on  ATS-6  on  day  243  of  1974. 

This  plot  displays  two  phenomena.  The  first  is  the 
ion  peaks,  which  are  clearly  presented  elsewhere  {Ref.  2}. 
These  ions  were  denoted  "spacecraft  generated  ions".  The 
second  is  an  apparently  analagous  phenomenon  in  the  electron 
spectrum.  These  observations  are  the  so-called  "Minnesota 
spots".  In  addition  to  the  electron  and  ion  spectra,  thi.^ 
figure  displays  the  detector  pitch  angle  in  the  top  section, 
which  is  related  to  the  rotation  angle. 

The  spots  are  seen  in  the  electron  spectrum  between 
100  and  300  eV  as  intense  white  points  that  appear  to  form 
triangle  shapes.  If  one  compares  the  periodicity  of  the 
electron  triangles  with  the  detector  look  angle  periodicity, 
a  direct  correspondence  is  evident.  As  the  top  of  the 
spacecraft,  where  the  detectors  are  located,  was  in  shadow 
at  this  time,  the  source  electrons  were  determined  to  be 
secondaries  from  the  surface  of  the  vehicle.  These 
electrons  had  been  emitted  from  a  surface  differentially 
charged   with   respect   to   the   detectors,    and   thus   the 
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Figure  17.  ATS-6,  Spectrogram,  Day  248  of  1974  {Ref.2} 
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charged  with  respect  to  the  detectors,  and  thus  the 
distribution  had  been  accelerated  by  a  potential  drop  <t> . 
The  source  of  the  flux  was  eventually  determined  to  be 
secondary  emissions  from  the  University  of  Minnesota 
detectors  also  located  on  the  satellite  measurements  box. 
{Ref.  2} 

3.   I  SEE 

Negative  charging  events  on  the  I3EE  satellites  were 
relatively  rare,  as  the  satellites  were  specifically 
designed  to  avoid  the  problem  of  spacecraft  charging.  ISEE- 
1  did,  however,  charge  to  significant  negative  levels  on  a 
few  occasions.  During  one  such  event,  analysis  of  the  ion 
mass  composition  revealed  some  interesting  features.  Figure 
18  is  a  presentation  of  the  ISEE-1  mass  spectrometer  reading 
from  March  17,  1973.  The  vertical  axis  is  a  count  rate,  and 
the  horizontal  axis  is  the  applicable  mass  channel.  The 
mass  channel  is  related  to  the  atomic  mass  by. 


(2-1) 


AMU  = 


3615 


Mass  Channel^ -^"^^ 


Thus,  mass  channel  54  corresponds  to  H"^,  15  to  N*,  and  14  to 
0*.  {Ref.  10} 

At  0533  UT  (Universal  Time),  the  instrument 
responded  with  a  broad  peak  in  the  high  Atomic  Mass  Unit 
(AMU=30-100)  channels.  This  peak  had  a  maximum  in  channel 
10  corresponding  to  32  AMU.  There  was  also  a  peak  at 
channel   54,   H*.    The  normal  value  of  the   ambient   atomic 
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Figure  18.  ISEE-1,  Mass  Spectrometer,  March  17,  1978 
{Ref.  10} 
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oxygen  is  indicated  by  the  solid  black  line  under  the  peak. 
The  data  from  0721  UT  are  quite  similar,  except  that  at  this 
point  the  maximum  has  shifted  to  channel  11,  28  AMU.  Again, 
this  is  larger  than  the  indicated  ambient  oxygen.  At  the 
time  these  observations  were  first  made  it  was  concluded 
that  these  molecules  were  probably  hydrocarbons  and  other 
outgassing  products.  In  light  of  spacecraft  surface 
sputtering,  however,  one  may  postulate  a  different  flux 
source.  That  is,  these  particles  are  ions  emitted  from  the 
spacecraft  surface.  Since  the  surface  is  primarily  silicon 
dioxide  (SiO-),  the  ejected  particles  would  largely  be  Si, 
0,  and  0-.  The  masses  of  these  are  28,  IG,  and  32  AMU 
respectively,  numbers  which  would  fit  nicely  into  the 
observed  data.  {Ref.l0} 

B.   P78-2  (SCATHA) 

Observations  of  spacecraft  generated  ions  are  presented 
for  five  days.  The  first  two  are  for  eclipse  charging 
events  in  1979,  the  second  pair  are  for  ion  gun  induced 
charging,   and  the  last  is  for  eclipse  charging  in  1931. 

1.   Day  83  of  1979 

This  first  eclipse  event  illustrates  the  triangle 
peak  and  its  associated  angular  dependence.  Figure  19  is 
the  grey  scale  spectrogram  for  a  one  hour  period  of  day  S3 
of  1979.  This  is  a  composite  of  the  electrons  from  the  Hi 
detector,  and  the  ions  from  all  three  detectors.  Each  of 
the   ion  displays  has  zero  energy  at  the  top   and   increases 


energy  downward.  The  electron  presentation  is  reversed, 
with  low  energies  at  the  bottom,  increasing  upward.  Between 
1740  and  1820,  the  triangle  structures  can  be  seen  in  the  Lo 
detector  data.  From  the  Hi  detector  data,  it  can  be  seen 
that  the  spacecraft  charges  negatively  to  about  200  v  at 
1750,  decreases  to  100  V  at  1815,  and  Increases  rapidly  to 
300  V  at  1825.  Multiple  shadow  peaks  are  visible  during 
this  period  in  both  the  Hi  and  Fix  detector  data. 

Figure  20  is  a  plot  of  the  log  ion  distribution 
function  (phase  space  density)  versus  energy  for  the  Hi 
detector.  The  spacecraft  is  charged  to  -250  V,  in  eclipse, 
as  indicated  by  the  peak  in  phase  space  density  (f)  at  250 
eV.  The  shadow  peak  is  the  secondary  peak  at  78  eV.  At 
this  time,  the  detector  is  parked  at  92  degrees, 
approximately    parallel   to   the   spacecraft    spin    axis. 

Figure  21  is  a  diagram  of  the  spacecraft  potential, 
triangle  peak  energy,  and  Lo  detector  look  angle  versus 
time.  The  energy  of  the  peak  varies  directly  with  the  look 
angle,  with  the  energy  minimum  occurring  at  the  angle 
maximum,  that  is,  looking  down  and  away  from  the  spacecraft 
(where  it  approaches  closest  to  the  spacecraft  surface). 
The  shadow  peak  energy  maintains  a  relatively  constant  ratio 
with  satellite  potential.  In  this  case,  Ep/Et  is  between  S 
and  7,  where  Ep  is  the  satellite  potential,  and  Et  is  the 
triangle  peak  minimum  energy.    This  is  significant,   a^   we 
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Figure    19.    SCATHA,    Spectrogram,    Day    83    of    1979 
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Figure  20.  Day  83  of  1979,  Ion  Distribution  Function 
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Figure  21.  Day  83  of  1979,  Vehicle  Potential,  Triangle 
Peak  Energy,  and  Look  Angle  vs.  Time 
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would  expect  a  differentially  charged  surface  to  maintain   a 
constant   ratio   over   time  with   the   satellite   mainframe. 
2.   Day  86  of  1979 

The  second  eclipse  charging  event,  day  36,  is 
similar  to  that  of  day  83.  Figure  22  is  a  grey  scale 
spectrogram  of  the  ion  data  for  a  one  hour  period  on  day  36 
of  1979.  The  triangle  structures  are  again  visible  from 
1645  to  1740  in  the  Lo  detector  display.  From  the  Hi  data 
we  see  that  charging  begins  at  approximately  1640,  and  the 
potential  varies  widely  throughout  the  period.  A 
fluctuation  in  the  satellite  potential  due  to  a  change  in 
the  electron  Alfven  boundary  is  seen  at  1705.  The 
spacecraft  achieves  a  maximum  potential  of  -500  V  at  1725. 
During  the  period,  a  shadow  peak  is  clearly  visible  at 
energies  lower  than  the  charging  peak  in  both  the  Hi  and  Fix 
detector  data . 

Figure  23  plots  the  ion  distribution  from  the  Hi 
detector.  The  satellite  potential  is  -340  V  as  indicated  by 
the  charging  peak.  There  is  a  well  defined  shadow  peak  at 
100  eV,  In  addition,  there  are  various  other  sm.all  peaks 
below  the  charging  peak  that  may  indicate  other  shadow 
peaks.  Frequently,  other  shadow  peaks  appear  in  the 
spectrograms,  maintaining  a  constant  separation  ratio 
between  each  other,  and  varying  as  the  main  peak  varier;. 
The  detector  is  again  parked  parallel  to  the  spacecraft  spin 
ax  is  . 
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Figure    22.    SCATHA,    Spectrogram,    Day    86    of    1979 
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The  strong  dependence  between  the  triangle  peak 
minimum  energy  and  spacecraft  potential  is  illustrated  by 
Figure  24  for  this  more  variable  day.  It  is  clear  that  the 
minimum  energy  of  the  triangle  peak  increases  with  the 
increasingly  negative  vehicle  potential.  The  ratio,  Ep/Et, 
is  maintained  between  6  and  G.5.  Again,  the  energy  of  the 
shadow  peak  varies  directly  with  the  detector  angle, 
achieving  minimum  energy  at  the  detectors  maximum  angular 
extent.  The  shadow  peaks  usually  lose  their  definition 
below  110  degrees. 

These  two  examples  illustrate  typical  observations 
during  eclipse  passages  for  which  negative  charging 
occurred.  By  comparison  with  ATS-6,  we  would  expect 
similar  observations  during  daylight  negative  charging 
events.  ATS-6,  however,  was  stabilized  and  had  its  large 
solar  array  to  provide  a  shadowed  area  which  would  simulate 
an  eclipse.  3CATHA  had  no  such  array,  and  in  addition  was 
spinning  so  that  no  area  was  shadowed  long  enough  to  charge 
to  a  high  level.  Observations  of  daylight  charging  showed 
that  the  spacecraft  potential  was  highly  spin  modulated, 
making  it  difficult  to  use  these  data  for  this  project. 
However,  induced  charging  events  with  the  ion  gun  provide  an 
additional  range  of  daylight  data.  Two  examples  of  such 
experiments  are  presented. 
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Figure  24.  Day  86  of  1979,  Vehicle  Potential,  Triangle 
Peak  Energy,  and  Look  Angle  vs.  Time 
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3.   Day  94  of  1979 

The  first  example  of  ion  gun  induced  charging  is  day 
94  of  1979.  The  ion  gun  was  in  operation  at  1  kV  and  1  mA . 
Figure  25  is  the  grey  scale  spectrogram  for  a  1  hour  period 
from  that  day.  The  data  presented  are  from  two  ion 
detectors,  with  the  Lo  detector  in  the  upper  diagram,  and 
the  Fix  in  the  lower.  The  ion  gun  is  turned  on  at 
approximately  1405  and  the  spacecraft  charges  promptly  to  a 
final  value  of  -G0  V.  Shadow  peaks  occur  at  11  and  20  eV  in 
the  Lo  detector  between  1420  and  1430,  and  between  1440  and 
1455. 

A  curious  feature  is  the  appearance  of  twin  shadow 
peaks  in  the  Lo  data  at  1440,  with  bands  between  in  which  no 
ions  were  detected.  When  this  phenomena  began,  the  detector 
had  just  been  parked  at  an  angle  of  20  degrees  to  the 
spacecraft  spin  axis,  away  from  the  satellite  body.  Txhe 
detector  had  previously  been  scanning.  When  the  data 
printouts  are  examined,  three  bands  clearly  emerge;  one  at 
11  eV,  another  at  22  eV,  and  the  last  at  30  eV.  We  may 
infer  that  certain  energies,  or  trajectories,  do  not  reach 
the  detector  at  that  angle.  These  events  indicate  something 
about  the  location  of  the  emitting  surfaces.  From  the 
striations  in  the  Lo  data,  we  would  expect  non-emitting  or 
non-d i f f erent ially  charged  surfaces  to  be  sandwiched  between 
emitting  regions,  thus  giving  the  dark  regions  where  no  ions 
are  detected.    Importantly,  the  Hi  detector  is  viewing  in  a 
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Figure  25.  SCATHA,  Spectrogram,  Day  94  of  1979 
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similar  direction  at  this  time,  and  gives  a  similar 
response.  The  Fix  detector,  viewing  70  degrees  from  the  Lo 
and  Hi,  reports  a  completely  different  environment.  It 
does,  however,  show  a  similarity  in  that  a  large  energy  gap 
exists.  This  indicates  mixed  differentially  and  non- 
different  ial  ly  charged  surfaces  as  above. 
4.   Day  200  of  1979 

Our  fourth  example  is  during  an  ion  gun  operation  at 
1  kV,  1  mA .  Figure  25  is  a  grey  scale  spectrogram  for  a  two 
hour  timeframe.  The  data  are  taken  from  the  Hi  detector. 
In  this  figure,  the  spacecraft  is  charged  to  approximately  - 
700  V,  except  during  the  two  five  minute  periods  when  the 
ion  gun  was  in  "trickle"  mode.  Trickle  mode  is  defined  as 
the  gun  in  operation  with  the  accelerating  voltage  off,  but 
current  still  applied  to  the  diode.  Thus,  the  generated 
ions  leak  out  at  low  energies.  Triangle  peaks  are  clearly 
visible  in  the  ion  display.  This  day  is  one  of  the  few  in 
which  these  peaks  become  evident  in  the  Hi  detector.  The 
ratio,  Ep/Et,  for  these  figures  is  in  the  9-10  range, 
indicating  that  the  Hi  detector  is  looking  at  an  emission 
point  different  from  that  seen  by  the  Lo  detector  on  days  83 
and  86  (or  different  processes  are  at  work). 

Figure  27  shows  the  distribution  functions  when  the 
Lo  and  Fixed  detectors  are  looking  approximately  parallel, 
radially  away  from  the  spacecraft.  The  Hi  detector  is  near 
its   maximum  rotational  angle.    The  satellite  potential   is 
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Figure  26.  SCATHA,  Spectrogram,  Day  200  of  1979 
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Figure  27.  Day  200  of  1979,  Ion  Distribution  Functions 
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approximately  -700  V.  The  Hi  and  Fix  detectors  exhibit  a 
broad  array  of  low  energy  ions,  but  such  a  spectrum  is  not 
evident  in  the  Lo  detector.  The  Lo  detector  does  show  a 
sharply  defined  shadow  peak.  The  data  have  a  similar  form 
throughout  this  period,  with  subtle  variations  in  the 
amplitude  and  energy  of  the  diffuse  spectrum,  particularly 
with  angle . 

5.   Day  92  of  1981 

Figure  28  is  the  grey  scale  spectrogram  from  a  1981 
eclipse  charging  event.  The  figure  displays  4  hours  of  data 
from  the  Hi  detector  for  both  ions  and  electrons.  The 
eclipse  charging  event  is  visible  in  the  illustration 
between  0810  and  0945  UT,  with  potentials  varying  from  -100 
to  -400  V.  During  this  period,  a  definite  shadow  peak  is 
also  evident,  in  the  -60  to  -150  V  region. 

Figure  29  is  a  plot  of  the  ion  distribution 
functions  for  day  92.  The  log  distribution  function  has  a 
charging  peak  at  0820  at  90  eV  with  a  value  of  6.07.  This 
value  is  comparable  to  those  on  days  83  and  86  of  1979. 
Thus  the  detected  countrates  are  as  high  as  those 
encountered  early  in  the  mission.  This  indicates  that  the 
source  of  the  detected  particles  is  not  dependent  on 
satellite  lifetime,  and  hence  that  outgassing  is  not  a  major 
contributor  to  the  observed  satellite  generated  ion 
phenomena . 
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Figure  28.  SCATHA,  Spectrogram,  Day  92  of  1981 
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Figure    29.    Day   92    of    1981,    Ion    Distribution   Function 
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Ill  .   DISCUSSION 

A.   POSSIBLE  ION  SOURCES 

The  ion  source  must  be  able  to  explain  three  different 
phenomena,  the  triangle  peaks,  the  shadow  or  mirror  peak, 
and  the  diffuse  spectrum.  The  prime  feature  that  the 
observations  have  in  common  is  the  appearance  of  the  subject 
ions  at  energies  below  that  provided  by  the  spacecraft 
potential.  Thus,  they  must  be  generated  in  regions  near  the 
spacecraft  where  they  do  not  experience  the  full  plasma- 
satellite  potential.  There  are  two  prime  candidates  for  the 
source,  outgassing  and  sputtering.  Outgassing  is  the 
emission  from  the  satellite  of  contaminant  molecules  and 
atoms  with  energies  near  the  temperature  of  the  satellite. 
These  are  primarily  trapped  hydrocarbons  and  other 
atmospheric  molecules.  The  outgassing  rate  may  be  increased 
by  the  inclusion  of  efflux  from  the  attitude  control  jets, 
and  the  ion  beam  system.  Sputtering  is  the  emission  of 
charged  and  neutral  atoms  from  the  satellite  surface  due  to 
ion  bombardment,  and  its  theory  has  been  described 
previously.  A  qualitative  discussion  of  both  possibilities 
follows,  followed  by  a  numerical  calculation  of  the  number 
of  ions  that  would  reach  the  detectors  due  to  Sigmund- 
Thompson  sputtering  from  differentially  charged  surfaces. 
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1 .   Outqassinq 

As  the  contaminant  particles  leave  the  spacecraft 
they  are  subject  to  ionization  from  photons,  and  collisions 
with  ambient  ions  and  electrons.  As  the  molecules  are 
ionized  at  various  distances  from  the  spacecraft  they  will 
return  in  a  broad  range  of  energies.  This  is  a  possible 
explanation  of  the  observed  diffuse  spectrum.  The  question 
is  then,  how  many  particles  outgas,  and  what  percentage  of 
these  are  ionized?  The  problem  of  outgassing  and  the 
outgassing  yield  has  been  studied  previously  {Refs.  32,33}. 
The  dominant  mechanism  for  the  ionization  of  the  outgassing 
molecules  has  been  determined  in  sunlight  to  be 
photoionizat ion,  and  in  eclipse  to  be  collision  with 
electrons.  The  electron  ionization  rate  is  approximately 
two  orders  of  magnitude  less  than  the  photoionizat ion  rate. 
In  eclipse,  however,  the  photoionization  rate  goes  to  zero 
and  the  electron  collisions  become  dominant.  In  all 
probability,  the  electron  ionization  rate  will  increase 
slightly  in  eclipse,  as  the  electron  density  will  increase 
due  to  space  charge  effects.  The  ionization  rate,  however, 
will  not  exceed  that  of  the  photoionization.   {Refs.  2,32} 

The  distance  within  which  the  ions  may  be  returned 
to  the  spacecraft  will  be  the  width  of  the  surrounding 
plasma  sheath.  The  sheath  distance  at  geosynchronous  orbit 
is  approximately  100  m.  With  this  information  and  the 
ionization  rates  above,   it  may  be  shown  that,  at  a  maximum. 
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approximately  .02%  of  outgassing  silicon  atoms  will  be 
ionized.  In  a  central  force  field,  with  conservation  of 
angular  momentum  and  energy,  calculations  reveal  that  all 
the  particles  ionized  in  the  sheath  will  be  returned  to  the 
spacecraft.  The  final  result  then,  is  that  .02%  of  the 
outgassing  molecules  vill  be  reattracted  to  the  satellite. 
{Ref.  32} 

When  the  spacecraft  is  first  launched,  the  ionized 
outgassing  flux  may  be  sufficient  to  provide  a  significant 
number  of  particles  at  the  detectors.  Outgassing,  however, 
is  an  extremely  time  dependent  phenomenon.  Once  the 
particles  leave  there  is  no  way  to  replenish  them,  so  the 
outgassing  rate  must  decrease  with  time.  Various 
experimental  data  have  placed  the  outgassing  as  proportional 
to  t~",  where  n  varies  from  0.8  to  1.6,  or  as  proportional 
to  exp( -t/41 . 6 ) ,  with  t  in  days  {Refs.  32,33}.  The  longest 
e-folding  time  then  would  be  41.6  days.  Even  this  would 
lead  to  a  decrease  in  the  outgassing  flux  by  three  orders  of 
magnitude  over  the  period  of  a  year.  As  noted  in  day  92  of 
1981  above,  however,  the  shadow  peaks  have  not  diminished 
over  time  in  any  significant  way.  This  leads  to  the 
conclusion  that  the  outgassing  flux  may  be  a  contributor, 
but  is  not  the  prime  source  of  the  ions. 

What  effect  do  the  thrusters  and  ion  gun  have  on 
this  conclusion?  An  experiment  was  performed  aboard  the 
SCATHA   spacecraft  to  measure  the  contamination   (impurities 
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mixed  with  the  hydrazine  fuel)  on  the  satellite  due  to  the 
thruster  operation.  Two  sensors  were  placed  on  the 
satellite,  one  in  the  bellyband,  and  one  on  the  forward  end. 
The  sensors  were  capable  of  detecting  5  ng/cm^  particle 
deposition  and  a  current  as  low  as  10"^^  A.  The  thrusters 
were  located  on  the  aft  end  of  the  vehicle,  and  were  fired 
once  per  week  during  the  bulk  of  the  flight  time  to  correct 
for  precession.  The  results  of  the  experiment  indicated 
that  no  measurable  flux  of  thruster  contaminants  was 
returned  to  the  spacecraft.  {Ref.34} 

The  ion  gun  is  also  an  unlikely  source  of 
contaminants  to  the  detectors.  As  noted  above,  it  is 
located  on  the  aft  end  of  the  spacecraft  and  directed  away 
from  it.  During  most  normal  modes  of  operation,  the  energy 
and  angular  momentum  of  the  exiting  particles  would  prohibit 
them  from  reaching  the  detectors.  There  is  evidence  that 
some  ions  are  returned  to  the  detectors,  but  they  are 
observed  at  energies  greater  than  the  satellite  potential 
{Ref.  35}.  There  is  a  possibility  that  ions  emitted  during 
trickle  mode,  at  low  energies,  may  be  reattracted,  and 
contribute  to  the  flux  observed  at  the  detectors.  This  is 
probably  only  a  small  intermittent  component,  as  the 
operation  of  the  ion  gun  in  the  above  mode  was  not 
continuous  or  frequent.  Thus,  the  outgassing  flux  is 
primarily  composed  of  the  emitted  contaminants. 
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With  just  the  emissions  of  the  contaminant  flux,  it 
is  still  numerically  conceivable  that  the  outgassing  flux 
could  produce  sufficient  particles  to  account  for  the 
observations.  However,  it  is  certain  that  the  outgassing 
flux  does  decrease  over  time,  and  this  is  not  observed  in 
the  data.  Additionally,  it  is  hard  to  conceive  that  the 
satellite  potenials  are  such  that  the  randomly  ionized  flux 
is  focussed  into  a  beam  one  energy  channel  vide  at  the 
detector.  Moreover,  it  is  unreasonable  that  the  focussing 
should  be  such  that  a  small  change  in  the  detector  look 
angle  yields  such  a  definite  change  in  the  observed  energy, 
as  in  the  phenomena  of  the  triangle  peaks.  In  sum  then,  it 
does  not  appear  likely  that  ionization  of  outgas  products  is 
the  source  of  the  observed  flux. 
2 .   Sputtering 

When  this  work  began,  it  was  initially  assumed  that 
the  ion  source  was  outgassing.  As  the  work  progressed,  it 
became  apparent  that  there  were  serious  drawbacks  in  this 
assumption,  as  noted  above.  While  searching  for 
alternatives,  the  advantages  of  sputtering  as  a  source 
became  clear.  First,  besides  a  small  dose  dependence,  there 
is  no  reason  for  the  sputtering  yield  to  decay  over  time. 
This  attribute  was  very  important  in  light  of  the 
significant  shadow  peaks  observed  in  1981.  Second,  the 
emission  of  sputtered  particles  from  a  differentially 
charged   surface   would   lead   to   observed   beams   at    the 
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detectors.  As  the  particles  left  the  surface  with  their 
various  initial  angles  and  energies,  they  would  be 
accelerated  by  the  existing  electric  field  into  an  almost 
parabolic  orbit.  Thus,  for  each  energy,  only  particles 
following  specific  trajectories  would  reach  the  detector, 
limiting  observation  angles.  Since  the  detector  energy 
channels  have  a  resolution  of  20  %,  the  bandwidth  of  the 
channels  increases  with  energy.  Therefore,  a  peak  which 
would  appear  broad  when  distributed  over  the  lower  energy 
channels,  can  be  compressed  into  a  narrow  peak  in  the  high 
channels.  The  low  and  high  energy  channels  of  the  SCATHA 
detectors  are  listed  in  Table  8.  The  similarity  between  the 
observed  data  and  the  data  from  the  Minnesota  spots  also 
became  evident.  It  seemed  reasonable  that  if  such  a  process 
could  occur  for  secondary  electrons,  it  could  occur  for  what 
were  essentially  "secondary"  ions.  Additionally,  the 
sputtering  theory  could  explain  the  diffuse  spectrum  as  the 
receipt  of  particles  sputtered  from  a  non  differentially 
charged  surface.  Without  the  charging,  the  ions  are  not 
focussed  into  the  higher  channels,  and  can  appear  as  a  broad 
low  energy  distribution.  The  third  phenomena,  the  shadow 
peak,  can  be  described  by  two  mechanisms.  The  first  is  a 
motionless  detector,  not  scanning  in  angle,  and  observing 
one  trajectory,  i.e.,  one  energy.  Alternatively,  the 
detector  could  be  scanning  laterally  across  an  emitted  beam. 


71 


as  opposed  to  vertically,  and  thus  be  picking  up  the  same 
energy  at  each  angle.   Figure  30  illustrates  the  geometry. 

The  problem  with  sputtering  as  a  source,  was  that  it 
had  never  seemed  to  be  able  to  provide  an  adequate  flux  of 
ions.  However,  the  only  previous  mention  of  sputtering  had 
been  prior  to  the  advent  of  the  Sigmund  theory,  and  before 
it  was  acknowledged  that  the  magnetospher ic  regions 
contained  large  numbers  of  oxygen  ions  {Ref.l}.  Preliminary 
calculations  at  this  point  indicated  that  the  flux  due  to 
oxygen  sputtering  was  sufficient  to  account  for  the  observed 
phenomena.  It  then  remained  to  model  the  situation  and 
determine  more  accurately  the  flux  at  the  detectors  due  to 
sputter  ing . 

It  was  now  necessary  to  evaluate  the  physical 
characteristics  of  the  environmental  interaction  with  the 
spacecraft.  As  the  bulk  of  the  vehicle  surface  was  covered 
with  solar  panels,  the  most  likely  target  of  the  incident 
ions  would  be  amorphous  silicon  dioxide.  For  a  magnetically 
disturbed  time  period,  the  ion  flux  would  be  composed 
primarily  of  oxygen  and  hydrogen.  The  Sigmund  theory  has  a 
strong  dependence  on  the  mass  of  the  incident  ions,  and  it 
is  reasonable  to  expect  a  much  larger  sputtering  yield  for 
the  heavier  ion.  See,  for  example.  Figure  14  for  H*  and  0"^ 
on  amorphous  silicon. 

The  next  step  was  to  use  the  equations  (1-7)  and  (1- 
8)   to  calculate  the  yield.    One  roadblock  to  this  was   the 
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Detector 

Scamng  Vertically  Through  Energy 


Detector  Scanning  Across  Angle  Receiving 
Constant  Energy 


Figure  30.  Detector  Beam  Viewing  Geometry 
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inability  of  the  Sigmund  theory  to  adequately  treat 
sputtering  from  mult i -atomic  targets.  There  are  some  data 
points  available  for  impacting  molecular  oxygen  on  silicon 
and  silicon  dioxide^  and  they  are  listed  in  Table  9 
{Ref.36}.  If  the  data  from  Table  7  are  applied  to  Table  9, 
we  see  that  the  sputtering  yield  for  0"^  on  silicon  will 
range  from  .1  to  .18  at  10  keV.  Similarly,  the  yield  on 
silicon  dioxide  will  range  from  .23  to  .44  atoms  per  ion  at 
10  keV.  The  sputtering  yield  calculated  from  the  Sigmund 
theory  for  atomic  oxygen  on  silicon  is  too  high,  and  more 
accurately  approximates  the  yield  from  SiOz.  We  will, 
however,  take  the  worst  case  approach,  and  calculate  the 
yield  for  silicon.  To  do  this,  we  will  increase  the  binding 
energy,  to  give  us  values  more  in  line  with  those  from  the 
experiment.  The  two  yield  curves  are  displayed  in  Figure 
31.   The  new  binding  energy  is  26.0  eV. 

Given  that  this  new  curve  now  reasonably  describes 
the  sputtering  yield,  the  problem  is  to  calculate  how  many 
particles  emitted  from  a  surface  will  reach  a  differentially 
charged  detector  for  a  realistic  incident  flux.  The  total 
yield  of  sputtered  particles  will  be  described  by 


(3-1) 


Y  = 


r^  s(v,e)  V 


f(v)  d^v 


In  this  expression,  r  is  the  maximum  energy  a  sputtered 
particle  can  have,  that  is,  the  maximum  amount  of  energy 
transferred.    This   maximum   is  calculated  by 
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Figure    31.    Oxygen-Silicon   Yield    Curves 
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r   =   4    M:l    M2       Em, 
(Mx+Mz)^ 


and  is  equivalent  to  the  maximum  energy  transferred  in  a 
binary  collision.  The  function  f(v)  is  the  distribution 
function  of  the  incident  ions,  and  S(v,9)  is  the  sputtering 
y ield  c 

With  the  total  yield  in  hand,  ve  may  now  attempt  to 
calculate  the  countrate  at  the  detector.  The  total  number 
of  counts  the  detector  receives  may  be  described  by  the 
following  equation. 


(3-2) 


/N 


Counts=Pdt  PdAPf  (V)  (\?.n)d^v 


This  assumes  that  the  detector  is  differential  in  time, 
area,  solid  angle,  and  energy.  If  we  integrate  over  the 
area  increment  and  time,  we  get, 

(3-3)      Counts  =  T  (3A  pdii  pv^  f  (  v  )  (  v- ri )  dv 


rr 


Thus,  the  countrate  is 


(3-4) 


Countrate  =  Counts/Time=3Arrv^f(v)(^.  n)dvdil 


If 


To  continue,  we  must  assume  that  limits  of  integration  on  d£ 
are  not  a  function  of  energy,  and  that  the  value  of  the 
integral  may  be  given  by  a  simple  number,  characteristic  of 
the  detector.  This  is  not  actually  true,  as  the  limits  on 
integration  depend  strongly  on  energy,  since  f(v)  may  not 
fill  the  entire  5  by  7  degree  aperture  of  the  detector.   The 
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inaccuracy  of  this  assumption  will  surface  later.    With  the 
assumption,  ve  derive. 


(3-5 


Countrate=CR  =  (5A  62    fv^  f  (  v  )  (  v- h  )  dv 


r 


The  two  constants,  (3A  and  6Sl^  are  combined  into  the 
geometric  factor  G,  a  constant  for  the  detector.  For  our 
detectors,  G  =  1.6  *  10~=*  m=  ster.  Due  to  the  differential 
aperture  of  the  detector,  the  value  of  the  dot  product  v-  h 
is  approximately  v  (i.e.  cos  9~1). 

We  must  now  determine  the  distribution  function  of 
the  emitted  particles.  From  the  earlier  discussion  of  the 
Thompson  Theory,  it  is  evident  that  the  energy  distribution 
of  the  yield  flux  is  proportional  to  E/(E+Ua)^.  That  is,  we 
have  the  function 


(3-6)      dY  =  A    E 

dE      (E+Ub)^ 


The  flux  Y,  is  also  given  by 


(3-7) 


Y  =  B  p«  E  f(E)  dE 


J 


where  f(E)  is  the  distribution  function  as  a  function  of 
energy  of  the  emitted  particles,  and  B  is  a  constant  that 
involves  the  integral  over  the  solid  angle.  If  we  combine 
these  two  relationships,  it  is  evident  (disregarding 
constants)  that 


(3-8)      f(E)  =  1  dY  =     1 

E  dE    (E  +  Ub)^ 


77 


We  will  now  determine  the  appropriate  normalization 
constants.  From  the  definition  of  the  distribution 
f unct  ion. 


(3-9) 


-r. 


f(v)  d^v 


Continuing  the  calculation,  we  assume  that  the  yield  flux 
and  the  distribution  function  have  similar  angular 
dependencies.  Thus,  E ( v, 9 ) =f ( v ) cose .  If  we  break  down  d^v, 
and  substitute  for  energy,  we  get 


(3-10) 


n  =  APd0Psin( 


e  cose  depf(E)i2Ei^  ^^ 


1 


In  the  case  of  emission  from  a  surface,  *  will  range  from  0 
to  2x1,  and  9  will  vary  from  0  to  k/2.  Integrating  over 
these  limits,  the  equation  reduces  to 


(3-11)     n  =  ATt(2/m3)^ 


J  (E  +  1 


dE 


Ub)^ 

This  integral  will  extend  the  full  energy  range  of  the 
emitted  particles,  that  is,  from  0  energy  to  r,  the  maximum 
energy  transferred  in  the  collisions.  If  one  then 
substitutes  a  dummy  variable,  x^,  for  E  the  integral  has  a 
standard  form  {Ref.  37}.   Solving,  we  then  get 


n=ATt(2/m3)^  P'' 


Ja   (E  +  U. 


dE 


n=ATt(8/m3)^ 


J 


dx 


s  =  r' 


<a   (X^  +  Ub)3 
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n  =  ATi(8/m=)^  {(-x)/(4(x2  +  U0  )  ^  )  +  x/(8U0(x2  +  Ua  )  ^  ) 
+  (l/8(Ua^)^)  tan-=^  (x/(U0)^)  }-a 

n=An:(8/m^)^  {(r^/(8(U0  +  T))){l/U<^    -    2/(Ue,  +  D) 
+  (l/(8Ua^^^))  tan-^  (r/Ua)^} 
Defining  f3  as  the  large  constant  on  the  right  side,  ve  get 

n=ATt(8/m^)"^  0 
or 

(3-12)     A  =  n  m  =  ^2/T^38^ 

Using  this  normalization  constant,  the  distribution  function 

becomes 

(3-13)     f(E,9)=  (n  m^/2/T^08v*)  ^^^q  1 


(E  +  Ue,)^ 
If  we  then  calculate  the  yield  flux,  Y,  ve  obtain 


Y=(l/2m)^   (n/13) 


P*"      E 

Jo  (E  +  U, 


dE 


)  = 


(3-14)     Y=(l/2m)^  (n/3){-(r  +  Ua/2)/(r  +  U0 ) =  +  (I/2U0)} 


or 


Y=(l/2m)^  (n/3)  6 


(3-15)     n=(2m)^0Y/a 

Therefore,  substituting  in  the  distribution  function. 


(3-16)     f(E,e)  =  ml  Y     1 

2k    6     (E  +  Ub) 


cos9 
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We  now  have  an  expression  for  the  distribution 
function  of  the  emitted  particles.  This  expression  can  then 
be  inserted  with  one  modification  into  eqn .  (3-5)  to 
determine  the  countrate  at  the  detectors.  The  necessary 
modification  accounts  for  the  differential  charge  between 
the  emitting  surface  and  the  detectors.  For  the  particles 
to  be  accelerated  to  the  detector,  the  surface  must  be 
charged  positively  relative  to  the  detectors,  therefore  * 
will  be  negative.  Since  E  will  map  to  E  +■  q*  (Liouville's 
Theorem),  we  will  obtain 


(3-17)     f (E,9,0)=  ml  Y 


cos9 


2Tt  (5  (  (E-*)  +  Ub)  3 
Here,   we  have  assumed  that  the  ions  are  singly  charged,  and 
<t>      will   be   converted  to  eV  from   V   on   this   basis.    The 
countrate  expression  is  then. 


(3-18)     CR=G(2/m=) 


P  f(E,e, 


<t>)    E  dE 


(3-19)     CR=2GY  cose 


J 


dE 


To  find  the  countrate,  we  must  now  integrate  for  each  of 
the  detector  energy  channels  over  its  energy  range.  As  this 
resolution  is  20%,  the  upper  and  lower  limits  are  1 . lEc  and 
0.9Ec,  where  Ec=  is  the  central  energy  of  the  energy  channel. 
The  integral  in  eqn.  (3-19)  is  a  standard  form,  and  the 
countrate  may  finally  be  written  as. 
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(3-20)     CR  =  2GY  cosQ{-(2.  2E.^  ^    U=,    -    <i>)     +  (  1  .  8E^  +  Ua  -  0)  } 
71(5         2(1. 1E=  +  Ua  -  0)=   2(.9E=  +  Ua  -  0)  = 

If  the  energy  channel  limits  are  both  less  than  0,   then  the 

countrate  over  that  channel  is  zero.   If  the  lover  limit   is 

below  0,   then  the  integral  is  taken  between  the  upper  limit 

and  0 . 

B.   CALCULATION  OF  THE  SPUTTERING  FLUX  AT  THE  DETECTOR 

Using  the  Sigmund  Theory  developed  above,  and  the 
derived  expression  for  the  countrate  at  the  detector,  we 
developed  a  program  to  calculate  the  sputtering  yield, 
incident  and  emitted  flux,  and  the  reponse  of  the  detector 
to  changes  in  both  the  energy  channel  and  differential 
charge,  0.  Data  for  the  ambient  plasma  compositions  was 
available  from  the  Lockheed  Ion  Mass  Spectrometer  for  days 
83,  86,  and  200  of  1979.  The  values  utilized  are  listed  in 
Table  10. 

The  first  portion  of  the  program  was  the  calculation  of 
the  total  yield  or  emitted  flux.  This  value  is  calculated 
from  the  following  expression; 

(3-21)     Y  =  P-  V  f(v)  S(v,9)  d=v 


I 


In  this  relation,  v  f(v)  is  the  flux  of  the  ambient  ions  as 
a  function  of  velocity,  and  S(v,9)  is  the  sputtering  yield 
as  a  function  of  velocity  and  angle.  The  factor  Sr.(£)  in 
S(v,e)  exists  in  tabulated  form,  and  is  not  described  by  a 
simple  continuous  function.    If  the  points  are  fitted  to   a 


curve,   ve  arrive  at  the  equation; 

(3-22)     S„(£)=  .35874  -  (.17885)log£  -(. 14359 )( log£ ) = 

+( .06193) (loge) =  +  ( .04228) doge) * 
-( .00541) (log£ )=  -(. 00630) (logE)« 
-( .00092) ( loge)^ 

This   equation   and   the   available   tabulated   points    are 

displayed   in   Figure   32.    Since  Sri(£)   is   a   complicated 

function,   the   integral   was   evaluated   numerically.    The 

function  was  calculated  at  a  thousand  points  over  an   energy 

range   that   depended   on  the   temperature   of   the   ambient 

Maxwellian  distribution.    For  a  4000  eV   distribution,   the 

high  energy  was  50000  eV  and  the  low  energy  was  100  eV.   100 

eV  was  chosen  as  the  lower  limit,   as  most  sputtering  yields 

are   approximately   zero  by  that  point.    At   the   50000   eV 

point,  the  value  of  the  distribution  function  has  reduced  by 

a  factor  of  10"^"  from  the  peak  value. 

The  yield  obtained  in  this  method  was  then  multiplied  by 

.02,   as  approximately  2%-4%  of  the  emitted  flux  is   ionized 

{Ref.   29}.    This  ion  yield  is  substituted  in  the  countrate 

integral   to  obtain  the  detector  responses.    The  values   of 

the    integrals    were   varied   over   energy   channel    and 

differential  charge  to  determine  if  the  responses  were   very 

narrow   as  observed.    Figures  33   shows  the  values   of   the 

incident,   emitted,   and  ion  fluxes  for  oxygen  as  a  function 

of  energy  for  day  86  of  1979.    The  figures  for  days  83   and 

200    are   similar.   From   this  it  is  very   clear   that   the 

incident   flux   is  much  greater  than  the  emitted   ion   flux. 
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Figure    32.    Lindhard's    Universal    Function    {Ref.27} 
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Thus  there  is  a  net  accumulation  of  positive  charge  on  the 
spacecraft  due  to  the  ambient  ion  current. 

Figures  34  and  35  are  displays  of  the  countrate  as  a 
function  of  energy  channel  for  a  constant  (35  eV) 
differential  charge,  for  days  86  and  200  for  oxygen  induced 
emission.  We  see  that  the  countrate  is  high  in  3-4  of  the 
energy  channels.  The  data  indicates  that  it  in  fact  should 
be  concentrated  in  only  one  channel.  However,  the  above 
equations  assume  that  the  detector  surface  is  a  flat  plate, 
and  that  all  particles  reaching  that  plate  are  counted. 
This  is  in  fact  not  true.  The  detector  viewing  cone  is 
probably  not  completely  filled,  as  noted  previously,  and  we 
are  seeing  the  error  due  to  treating  the  factor  <3il  as  a 
constant  in  eqn .  (3-3).  This  factor  will  serve  to  reduce 
the  countrate  overall,  and  since  particles  at  the  edges  of 
the  viewing  cone  (with  different  energies  from  the  cone 
center)  will  be  reduced,  the  peak  will  most  likely  become 
narrower.  The  response  for  a  specific  channel  and  varying 
differential  charges  is  quite  similar.  That  is,  a  given 
channel  will  only  respond  for  a  certain  narrow  range  of 
differential  charge.  Similar  calculations  for  the  ambient 
hydrogen  population  on  the  same  days  indicate  that  it  can 
also  generate  a  significant  flux  at  the  detectors,  simply 
due  to  the  large  amount  present.  The  relative 
contributions  of  these  two  ambient  plasma  components  very 
widely,   depending   on  the  temperature  and  density  of   their 
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Figure  34.  Channel  Response  for  a  35  V  Differential  Charge 
0+  Incident,  Day  86,  1979 
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Figure  35.  Channel  Responses  for  a  35 
0+  Incident,  Day  200,  1979 
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respective  distributions.  In  general,  the  oxygen 
contribution  is  dominant. 

The  countrates  from  the  detectors  are  available  for  the 
above  three  days.  On  day  83  of  1981,  the  measured 
countrates  varied  from  29  -  800  counts  per  second,  with  the 
average  at  360.  For  day  86,  the  range  was  24  -  1080,  with 
the  average  at  300.  Day  200  was  much  lower,  between  25  and 
213,  with  an  average  value  of  120.  The  available  mass 
spectrometer  data  was  averaged  over  one  hour  periods,  so  we 
will  generally  calculate  only  one  representative  number  for 
the  charging  time  frame.  The  calculated  countrates  were 
380,  325,  and  400  for  days  83,  86,  and  200,  respectively. 
Considering  the  approximate  nature  of  the  calculations,  the 
theoretical  responses  compare  favorably  with  the  measured 
data,  with  the  exception  of  day  200.  However,  it  has 
already  been  noted  that  this  data  from  the  Hi  detector 
indicated  a  sputtering  source  different  from  those  seen  on 
the  Lo  detector.  This  further  supports  that  observation,  as 
the  sputtering  yield  for  this  event  is  clearly  much  lower. 

As  the  assumptions  in  this  calculation  tend  to  minimize 
the  number  of  detected  particles,  it  seems  reasonable  to 
conclude  that  sputtering  from  a  differentially  charged 
surface  is  sufficient  to  account  for  the  narrow  triangle 
peaks.  These  discrete  responses  will  serve  to  explain  the 
triangle  and  shadow  peaks.  This  leaves  the  diffuse 
background.   We  require  a  broad  response,  with  the  particles 


scattered  equally  across  the  energy  channels.  Figure  36 
gives  the  log  distribution  function  as  a  function  of  energy 
for  an  emission  angle  of  0  degrees.  The  distribution 
function  has  its  maximum  in  the  low  energy  channels,  and 
numerically  compares  favorably  with  Figures  21  and  24.  The 
broad  spectrum  may  then  be  described  as  the  returning 
emissions  from  non-d i f f erent ially  charged  surfaces. 

We  must  now  explain  the  shadow  peaks.  These  appear  to  be 
a  combination  of  two  occurrences,  acceleration  from  a 
differentially  charged  surface  to  a  motionless  detector,  or 
acceleration  from  a  differentially  charged  surface  to  a 
detector  scanning  laterally  across  the  emitted  beam,  as 
opposed  to  vertically.  See  Figure  31  for  a  view  of  the 
geometry.  Since  this  is  then  a  variation  of  the  triangle 
peaks,  the  source  and  results  would  be  the  same. 

The  trajectories  to  the  detectors  can  be  recreated  by 
the  NASCAP  program.  Figure  37  displays  the  particle 
trajectories  for  the  observed  angles  and  energies  at  the  Lo 
detector  on  day  86.  From  top  to  bottom,  the  lines  are  as 
follows;  a)84  eV,  145  degrees,  or  12.5  degrees  above  the 
spacecraft  plane,  b)  70  eV,  152  degrees,  c)  60  eV,  175 
degrees,  d)  and  53  eV,  188  degrees.  One  can  see  that  the 
particles  track  back  to  the  spacecraft  surface  to  different 
points.  These  points  are  all  located  in  the  solar  array 
portion  of  the  spacecraft  surface.  These  results  are  as 
expected.    The  solar  array  has  charged  differentially,   and 
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Figure  36.  Channel  Response  for  a  Zero  V  Differential 
Charge 
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Figure  37.  NASCAP  Trajectories,  Day  86  of  1979 
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the  sputtered  particles  follow  specific  trajectories  in  the 
electric  fields.  The  limited  viev  cone  of  the  detector  only- 
picks  up  specific  energies. 

NASCAP  was  also  run  for  the  potentials  observed  on  day 
94.  The  code  was  run  to  simulate  the  responses  of  the  Lo, 
Hi,  and  Fix  detectors  to  a  range  of  incoming  particle 
energies.  The  energies  ranged  from  2,0  to  90.1  eV.  As 
noted  above  in  the  observations,  the  spacecraft  was  charged 
to  -60  V,  and  peaks  were  noted  at  various  lower  energies  on 
the  detectors.  The  Hi  detector  was  parked  at  67  degrees, 
and  the  Lo  at  91  degrees.  Thus,  Hi  and  Lo  detectors  were 
sampling  approximately  the  same  region  of  space. 

Figures  38  and  39  depict  the  Lo  detector  trajectories  in 
two  different  planes.  The  first  shows  the  top  of  the 
spacecraft  surface.  Particles  at  energies  less  than  43  eV 
are  returned  to  the  spacecraft.  The  trajectories 
corresponding  to  observed  particles  are  marked  with  an 
asterisk.  The  second  figure  is  a  complementing  view  from 
the  left  side  of  the  first  figure.  These  two  figures 
indicate  that  the  actually  observed  particles  are  emitted 
from  the  satellite's  forward  end.  Additionally,  there  are 
gaps  from  which  no  particles  are  seen  at  the  detector.  An 
explanation  for  this  is  the  existence  of  patches  of 
different  nonconducting  material  on  the  forward  surface  for 
use  in  the  Satellite  Surface  Potential  Monitor  (SSPM).  When 
these  patches  (kapton,   silvered  teflon,   and  quartz)  charge 
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Figure  38.  NASCAP  Trajectory,  Day  94,  Lo  Detector 
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Figure  39.  NASCAP  Trajectory,  Day  94,  Lo  Detector 
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differentially,  the  particles  reaching  the  detector  from 
each  will  follow  different  trajectories,  and  have  different 
energies.  Additionally,  the  S-band  antenna  mast  on  the 
forward  end  is  nonconducting,   and  may  be  a  particle  source. 

Figures  40  and  41  are  similar  views  for  a  model  of  the 
Hi  detector  trajectories.  We  see  that  the  origin  of  the 
particles  is  also  on  the  upper  surface,  and  the  data  shows 
the  same  striations.  Again,  the  only  nonconducting  surfaces 
on  the  forward  end  are  the  SSPM  and  the  antenna.  The  detail 
on  NASCAP  is  not  sufficient  to  study  the  extremely  strong 
local  electric  fields  expected  in  this  region.  If  the 
satellite  potential  is  100  V,  then  the  field  due  to  the 
sheath  is  approximately  1  V/m.  In  the  forward  surface 
vicinity,  with  a  differential  potential  of  30  V,  the  fields 
will  be  on  the  order  of  30-60  V/m,  and  will  significantly 
distort  the  local  trajectories  of  particles.  Thus,  the 
NASCAP  trajectory  plotting  routine  is  somewhat  unreliable  in 
this  region.  Looking  at  both  figures,  however,  the  origins 
for  the  two  detectors  are  reasonably  close,  and  allowing  for 
the  local  fields,  may  possibly  come  from  the  same  surface. 
As  a  working  hypothesis,  we  will  assume  that  the  SSPM  is  the 
source  of  the  emitted  particles. 

Figures  42  and  43  model  the  fixed  detector.  The  views 
are  similar  to  those  above.  We  see  that  the  particles 
viewed  by  the  fix  detector  come  from  completely  different 
regions   of  the  spacecraft.    All  the  obviously   reattracted 
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Figure  40.  NASCAP  Trajectory,  Day  94,  Hi  Detector 
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Figure  41.  NASCAP  Trajectory,  Day  94,  Hi  Detector 
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Figure  42.  NASCAP  Trajectory,  Day  94,  Fix  Detector 
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Figure  43.  NASCAP  Trajectory,  Day  94,  Fix  Detector 
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particles  have  energies  of  less  than  43  eV,  as  did  the  other 
two  detectors.  We  can  see  that  the  trajectories  of  the 
particles  skip  over  the  bellyband  region.  This  is 
reasonable,  as  the  bellyband  would  be  expected  to  be  at  the 
satellite  mainframe  potential,  and  thus  particles  emitted 
would  not  be  accelerated  by  a  differential  charge.  The 
trajectory  calculations  for  the  Fix  detector  also  support 
the  theory  of  sputtering  from  differentially  charged 
sur  faces . 
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IV.  CONCLUSIONS 

It  is  reasonable  to  conclude  that  sputtering  is  a  likely 
source  of  the  observed  phenomena.  It  can  satisfactorily 
explain  the  flux  and  structure  of  the  triangle  and  shadow 
peaks.  It  also  can  account  for  the  observed  broad  spectrum 
of  ions.  These  fluxes  are  large  enough  to  be  measured  by 
the  detectors,  but  are  they  large  enough  to  cause 
substantial  damage  to  the  spacecraft  surface? 

For  this  data  set,  the  sputtering  yield  was  less  than 
the  incident  flux,  so  that  the  net  result  was  particle 
accretion  on  the  surface.  Coincidentally,  the  data  in  this 
thesis  were  taken  from  relatively  low  energy  days.  On 
higher  energy  days,  the  sputtering  yield  increases,  and  the 
net  result  is  surface  depletion.  However,  the  surface 
damage  caused  by  this  process  appears  to  be  minor.  Assuming 
a  higher  energy  day,  with  emission  on  the  order  of  twice  the 
incident  flux,  one  atomic  layer  will  be  removed  from  a 
silicon  dioxide  surface  every  8  years.  Since  the  solar 
cells  are  covered  by  the  glass  shields,  they  will  be 
impervious  to  the  surface  sputtering. 

Another  facet  of  this  is  the  damage  to  surface 
conducting  coatings  on  spacecraft.  These  coatings  are  only 
several  atoms  thick  and  may  be  expected  to  sputter  from  the 
surface  within  several  months  of  the  satellite  launch.   This 


101 


may  account  for  the  observed  charging  on  the  ISEE-1 
satellite . 

Since  the  energy  spectrum  of  the  sputtered  particles 
dies  off  quite  rapidly,  and  the  observed  differential 
charges  are  usually  not  large,  it  is  not  likely  that  any 
great  contamination  of  the  charging  peak  results.  The 
sputtering  in  general  will  probably  be  uniform  over  a  given 
material.  Surfaces  will  be  worn  down  over  large  areas  as 
opposed  to  experiencing  significant  localized  damage. 
Therefore,  it  appears  that  no  short  term  damage  will  result 
from  sputtering  on  spacecraft. 

Further  work  is  required  to  accurately  calculate  the 
sputtering  yields  from  various  surfaces  to  prove  or  disprove 
these  conclusions.  At  this  point,  given  the  current  state 
of  the  sputtering  field,  it  is  not  possible  to  state  with 
any  certainty  whether  extensive  damage  will  result  from 
ambient  ion  bombardment.  There  appears  to  be  no  work  in 
progress  anywhere  to  determine  experimentally  the  sputtering 
yields  for  most  ion-target  combinations.  Since  the  Sigmund 
Theory  is  so  limited,  these  results  are  necessary  for  any 
actual  calculations  of  the  effect  of  the  ambient  heavy  ion 
bombardment   on  space  vehicles. 

This  thesis  demonstrates  that  sputtering  is  the  probable 
cause  of  the  observed  low  energy  ion  fluxes.  Additional 
work  is  required  to  refine  the  NASCAP  trajectory  information 
and   couple  it  with  the  theory  developed  here  to   model   the 
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detector  response  accurately.  It  is  recommended  that  this 
work  be  continued  when  and  if  further  applicable  sputtering 
data  become  available. 
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APPENDIX 

TABLE  1 

MAGNETOSPHERIC  IONS  AND  THEIR  SOURCES  {Ref.  7} 

Origin Source  Ions   ._ 

Sun            Solar  Wind  U*,Ue**,0^* 

Earth          Ionosphere  H*,He*,0* 

Earth          High  Altitude  H*, He*, He**, 0*' 

Thermal  Plasma,  o**,0^* 
Plasmasphere 
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TABLE  2 


STORM  TIME  RING  CURRENT  COMPOSITION  {Ref.  8} 


Date 

07  July 

26  Nov 

09  Mar 

30  May 

1977 

1977 

1978 

1978 

H*  Flux 

la-'cm-^s- 

■1 

1.26 

4.20 

13.0 

6.0 

Dens  i ty 

10-=cm-3 

H* 

11 

44 

143 

58 

He* 

.72 

4.3 

10.9 

6.6 

0* 

8.3 

91 

228 

107 

Mean  Ener 

gy 

keV 

H-^ 

8.3 

5.7 

5.1 

6.4 

He* 

4.5 

4.3 

5.4 

2.7 

0* 

4.9 

3.8 

5.3 

4.  4 
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TABLE  3 

STORM  TIME  MAGNETOSPHERE  COMPOSITION  {Ref.  8} 

Date      04  May  29  July    29  Jan  02  June 

1977  1977  1978  1978 

H*  Flux   3.3  2.4  10  2.8 
10^cm-^s-^ 

Density 
10-=cm-^ 

H*    26  20  10  27 

He**  1.5  .16  9.3  .37 

He*   1.8  .9  .3  .79 

0*    3.0  13  4.9  60 

Mean  Energy 
keV 

H*    9.6  8.5  4.2  6.8 

He**  20.3  15.9  8.6  14.9 

He*   4.8  3.9  6.4  4.8 

0*    4.5  4.7  5.4  4.1 
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TABLE  4 

VALUES  OF  THE  REDUCED  STOPPING  POWER  {Ref.27} 
Sr,(£) e Sn(£) 


0.120 

0.002 

0.403 

0.20 

0.154 

0.  004 

0.  405 

0.  40 

0.211 

0.01 

0.356 

1.00 

0.261 

0.02 

0.291 

2.00 

0.311 

0.04 

0.214 

4.00 

0.372 

0.10 

0.128 

10.0 

0.393 

0.15 

0.0813 

20.0 

0.0493 

40.0 
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TABLE  5 


CALCULATED  AND  MEASURED  SPUTTERING  YIELDS  FOR  VARIOUS 
MATERIALS  {Ref.  23} 

THE  NUMBERS  IN  PARENTHESES  ARE  THE  YIELDS  CALCULATED 
FROM  THE  SIGMUND  THEORY 


TARGET 

UB(eV) 

Z^ 

SDutterinq  Ra 

t io  ( ate 

Dms/ion ) 

MATERIAL 

Ne* 

Ar* 

Kr* 

Xe^ 

Pb 

2.01 

82 

3.6 

10.5 
(29.7) 

24.0 
(44.0) 

44.  5 
(74.6) 

Ag 

2.94 

47 

4.5 

10.8 
(12.7) 

23.5 

(20.1) 

36.2 

(24.0) 

Sn 

3.11 

50 

1.8 

4.3 
(12.5) 

8.5 
(19.5) 

11.8 
(24.9) 

Cu 

3.46 

29 

3.2 

6,8 

11.8 

19  ,0 

(3.7) 

(6.7) 

(11.9) 

(15.5) 

Au 

3.79 

79 

3.6 

10.2 
(16.1) 

24.5 
(23.9) 

39.0 
(27.8) 

Pd 

3.87 

46 

2.5 

5.3 
[9.24) 

10.5 
(14.1) 

14.4 
(18.1) 

Fe 

4.29 

26 

1.3 

2.3 
:5.15) 

4.0 
(8.92) 

4.9 
(11.7) 

Ni 

4.43 

28 

1.4 

3.5 
5.35) 

5.6 
(9.01) 

7.6 
(11.8) 

V 

5.3 

23 

0.3 

1.0 
3.85) 

1.7 
(6.21) 

1.9 
(8.9) 

Pt 

5.82 

78 

1.9 

5.3 
9.35) 

11.3 
(15.3) 

16.0 
(19.2) 

Mo 

6.82 

42 

0.6 

1.5 

5.14) 

2.7 
(7.70) 

3.8 
(10.1) 

Ta 

8.06 

73 

0.7 

1.6 
6.62) 

3.1 
(10.0) 

4.0 
(12.3) 

W 

8.70 

74 

1.0 

2.3 
5.07) 

4.7 
(9.31) 

6.4 
(11.9) 
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TABLE  6 


CALCULATED  SPUTTERING  YIELDS  USING  A 
MODIFIED  BINDING  ENERGY. 


TARGET 

SDutter  inq 

Ratio 

( Atoms/I  on ) 

MATERIAL 

1.7  Ua 

Ar-^ 

Kr^ 

Xe^ 

Pb 

3.42 

17.5 

25.9 

43  .9 

Ag 

5.00 

7.47 

11.8 

14  .1 

Sn 

5.29 

7.35 

11.5 

14.7 

Cu 

5.88 

3.94 

7.00 

9.12 

Au 

6.44 

9  .47 

14.1 

16.35 

Pd 

6.58 

5.44 

8.3 

10.7 

Fe 

7,29 

3.03 

5.25 

6.9 

Ni 

7.53 

3.15 

5.3 

6.94 

V 

9.01 

2.26 

3.55 

5.24 

Pt 

9.89 

5.50 

9.00 

11.24 

Mo 

11.59 

3.02 

4.53 

5.94 

Ta 

13.7 

3.89 

5.88 

7.24 

W 

14.79 

2.98 

5.48 

7.00 
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TABLE  7 

RATIO  OF  SPUTTERING  YIELD  PER  ATOM  FOR  MOLECULAR 
AND  ATOMIC  ION  BOMBARDMENT  {REF.  31} 

Projectiles  Si  Ag  Au 

CI-CI2  ...  1.09 

Se-Sez  1.15  1.44  1.44 

Te-Te^  1.30  1.67  2.15 
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TABLE  8 


SCATHA  PARTICLE  DETECTOR  ENERGY  CHANNELS  {Ref.  4} 


Lo 

Detector 

Step 

Energy ( eV) 

Step 

Energy ( eV ) 

0 

-0.334 

32 

26.26 

1 

-0.29 

33 

30.10 

2 

-0.23 

34 

34.57 

3 

-0.16 

35 

39  .70 

4 

-0.08 

36 

45.55 

5 

0.00 

37 

52.22 

6 

0.09 

38 

60.00 

7 

0.22 

39 

68.77 

8 

0.  34 

40 

78.55 

9 

0.50 

41 

90.10 

10 

0.67 

42 

103. 32 

11 

0.87 

43 

118.32 

12 

1.10 

44 

135.65 

13 

1.36 

45 

155.43 

14 

1.66 

46 

178.09 

15 

2.00 

47 

203.98 

16 

2.39 

48 

233.42 

17 

2.83 

49 

267.20 

18 

3.34 

50 

306.30 

19 

3.94 

51 

350.74 

20 

4.60 

52 

401.96 

21 

5.37 

53 

457.73 

22 

6.25 

54 

525.50 

23 

7.25 

55 

601.05 

24 

8.40 

56 

687.71 

25 

9.72 

57 

787.70 

26 

11.23 

58 

902.13 

27 

12.96 

59 

1033.23 

28 

14.96 

60 

1184.33 

29 

17.23 

61 

1355.42 

30 

19.83 

62 

1529.85 

31 

22.80 

63 

1777.60 
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TABLE  8 

CONTINUED 

Hi  Detector 

step 

Energy  { 

eV)          Step 

Energy  (eV) 

0 

-4,68 

32 

1205.49 

1 

-2.35 

33 

1383.30 

2 

0.32 

34 

1586.90 

3 

3.38 

35 

1820.02 

4 

6.88 

36 

2086.94 

5 

10.89 

37 

2392.57 

6 

15.49 

38 

2742.51 

7 

20.75 

39 

3143.19 

8 

25.77 

40 

3601.97 

9 

33.66 

41 

4127.28 

10 

41.56 

42 

4728.76 

11 

50.60 

43 

5417.45 

12 

60.95 

44 

6206.00 

13 

72.80 

45 

7108. 89 

14 

86.37 

46 

8142.70 

15 

101.91 

47 

9326.41 

16 

119.71 

48 

10681.76 

17 

140.08 

49 

12233.63 

18 

163.41 

50 

14010.53 

19 

190.12 

51 

16045.08 

20 

220.70 

52 

18374.64 

21 

255.72 

53 

21041.98 

22 

295.81 

54 

24096.09 

23 

341.72 

55 

27593.05 

24 

394.29 

56 

31597.06 

25 

454.45 

57 

36181.66 

26 

523.39 

58 

41431.02 

27 

602.30 

59 

47441. 54 

28 

692.65 

60 

54323.58 

29 

796.11 

61 

62203.53 

30 

914.56 

62 

71226.06 

31 

1050.19 

63 

81556.86 
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TABLE  9 


THE  SPUTTERING  YIELDS  OF  10  keV  ION  BEAMS  ON 
SILICON  AND  SILICON  DIOXIDE  {Ref.  36} 


Target  Incident  Ion  Beam 

Ar*  O2* 

Si  0.94  0.20 

SiOz  0.69  0.48 
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TABLE  10 


AMBIENT  PLASMA  VALUES  FOR  DAYS  83,  86,  AND  200  OF  1979 
DATA  FROM  THE  LOCKHEED  ION  MASS  SPECTROMETER. 


(m~^)      Energy(keV) 

8.852 

120.778 

33.992 

4.217 

5.889 

43.728 

3.662 

4.035 

5.271 
7.279 
0.420 
1.957 


Day 

Particle 

Number 

Dens  i 

83/ 

H  + 

3.845 

*  10» 

1979 

Hs^ 

2.487 

*  10^ 

He 

2.865 

*  10  = 

0* 

3.188 

*  10  = 

86/ 

H* 

2.606 

*  10» 

1979 

Ha-" 

2.850 

*  103 

He 

2.410 

*  103 

0^ 

1.635 

*  10» 

200/ 

H^ 

9.064 

*  10^ 

1979 

H2* 

4.209 

*  10* 

He 

2,169 

*  10  = 

0* 

4.093 

*  10  = 
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